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CHAPTER 1 



Introduction 



The problem of selective excitation pulse design in nuclear magnetic resonance (NMR) corre- 
sponds, in mathematics, to the inversion of a certain mapping, T, which we call the selective exci- 
tation transform. This transform, which is a kind of non-linear Fourier transform, maps complex- 
valued functions of time, called pulses, to unit 3-vector valued functions of frequency, called mag- 
netization profiles. The explicit definition of T is given in Section [LT] The theory of inverting T 
has been shown to coincide with the theory of inverse scattering for the Zakharov-Shabat (ZS) 2x2 
system (see Section IZ21 and for example |14|). Numerous authors have studied this inverse scatter- 
ing problem, for example Ablowitz et al. 1 1 ] , and the results have been applied to NMR pulse design 
(see for example |[l4l l3l). However, no stable and efficient algorithm has been given in the literature 
for generating the full space of solutions to the inverse problem. For this and other reasons, less 
exact methods for NMR pulse design, such as the Fourier transform method and the Shinnar-Le 
Roux (SLR) method, have been used in practice instead of the more exact and more flexible inverse 
scattering (1ST) method. In this thesis, we present the discrete inverse scattering transform (DIST) 
algorithm for efficiently solving the full inverse scattering problem relating to NMR pulse design. 

In this introductory chapter, we define the continuum and discrete selective excitation trans- 
forms and describe the problem of selective excitation pulse design. The theoretical results are 
summarized in Sections [T31 II. 61 and ll.7[ and the main algorithms are described in Section [L8l 

The second chapter is devoted to proving the main theorems and deriving the algorithms. We 
introduce the discrete scattering theory, which is completely analogous to the standard continuum 
theory. 

The third chapter describes how to apply the theory to practical NMR pulse design. 
See [4J and [5 J for a detailed mathematical introduction to NMR imaging and NMR pulse 
design. 



In this section we define the selective excitation transform, which maps a complex function of 
time to a unit 3-vector valued function of frequency. 

Let a; : M ^ C be a function of time (usually we think of u as smooth and supported on a finite 
interval), and suppose that for every frequency z G M, there is a solution M_(z; •) : R — > to the 
frequency dependent Bloch equation (without relaxation) 



1.1. The selective excitation transform 



(1.1.1) 



^ [Re Lo{t) 

—M_{z;t) = M^{z;t) X Imu;(t) 



normalized by 



(1.1.2) 



t— >— oo 



lim M_(z;t) 
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Let us show that such a solution is unique. If Mi and M2 are two solutions to (ll.l.ll) then 

Mi{z; t)'^M2{z; t) is independent of t because 



d 



-Miiz;tfM2iz;t) 



[Xb (z; t)Mi {z; t)f M2 {z; t) + Mi{z; tf [Xb (z; t)M2 (z; t)] 



Mi 



z;t) {XB{z;ty + XB{z;t)) M2{z;t) 



where 

z — ImLi;(t) 

XB{z;t)= -z Reuj{t) 

lmuj{t) -Reuj{t) 

Therefore, (I1.1.2I ) implies that M_(z; t)^M_(z; t) = 1. So M^{z; t) is a unit vector for all z and 
t. It is unique because, if M2{z;t) is any unit vector satisfying M-{z;t)'^ M2{z;t) = 1, then the 
Cauchy-Schwarz inequality implies that M_(z; t) = M2{z; t). 
Suppose that u decays sufficiently so that 

Re e*^* — Im e*^* 



M{z) :-- 



lim 



Im e*^* 




Re e 




izt 



M^{z;t) 



exists for all t. For example, this limit exists whenever lo is integrable (see S). We call to the 
pulse and M the resulting magnetization profile. The map a; M is called the selective excitation 
transform, and we write M = Tu>. In Section [T3] we plot several examples of pulses and their 
resulting magnetization profiles. 



1.2. The problem of selective excitation pulse design 

In selective excitation pulse design, we usually start with an ideal magnetization profile Mijeai : 



c 



This is a unit 3-vector valued function of frequency which is typically equal to 



[0 1 ] outside some finite interval. The problem is to find a pulse a; : M ^ C such that 
the resulting magnetization profile M = Tuj is a good approximation to Mjdeai- For practical 
applications, uj should have finite duration. That is, it should be supported in some finite interval 
[p — T, p] C M. The number T is called the duration of the pulse, and the value p is called the 
rephasing time. Depending on the application, it may be important to design a pulse with the 
shortest possible duration and minimal (perhaps zero) rephasing time (pulses with p = are called 
self refocused). It is also often necessary to limit the energy 



E.., 



\uj{t)Y'dt 



of the pulse, as well as the maximum amplitude. For many application it is important that the pulse 
behaves well under imperfect magnetic field conditions. That is, should not be too large. 
The standard example of an ideal magnetization profile is 



Midealiz) 





sin ^0 
cos 60 






if I z — zq I < Co 



if \z — Zq\ > Cq. 
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Figure 1.3.1. A constant pulse and the resulting magnetization: T = I, p = |. 
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Here, the magnetization is rotated 9o radians around the x-axis for frequencies in the interval | z — 
zo\ < cq. Outside of this frequency interval, the magnetization is kept at equilibrium. The angle is 
called the flip angle. See the next section and Chapter[3]for plots of pulses producing magnetization 
profiles that approximate such a profile. 



1.3. Examples 



When \uj{t)\dt is small, the resulting magnetization profile M = Tu very closely resem- 
bles the Fourier transform of uj (see for example |5|). This fact is illustrated in Figures [T.3.1l and 
11.3.21 Figure 11.3.11 shows a pulse which is constant over a finite interval. Notice that at zero fre- 
quency, the magnetization is rotated ^ radians around the x-axis, and away from zero the transverse 
magnetization resembles a sine function. Figures [T.3.2l and ri. 3. 3l show truncated sine pulses. Such 
pulses are often used in practice to produce a magnetization profile which is approximately con- 
stant within some finite interval, and approximately in equilibrium outside the interval. Notice that 
as the flip angle increases from 90° to 120°, the resulting magnetization profile does a worse job 
approximating the ideal. 

To obtain a more accurate magnetization profile, one should invert the selective excitation trans- 
form directly rather than use the Fourier transform approximation. Figures ll.3.4[ ll.3.5[ and 11.3.61 
show three very different pulses which produce magnetization profiles which accurately approx- 
imate a 120° profile. The fact that the inversion of the selective excitation transform is highly 
non-unique reflects the nonlinearity of the map uj i— > Tlo. 



7 




8 



Figure 1.3.4. A 120° minimum energy pulse and the resulting magnetization. 
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Figure 1.3.5. A 120° self refocused pulse. 
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Figure 1.3.6. A 120° self refocused pulse. 
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1.4. The discrete selective excitation transform (hard pulse approximation) 

For practical purposes it is useful to consider pulses of the form 



j=-oo 



For example, this is the type of pulse designed by the SLR algorithm (see [12]). Such a pulse is 
called a hard pulse. If ujj vanishes for all j > p ^ Ia, then Q, is said to have p rephasing time steps. 
In practical applications, a softened version, e.g.. 



(1-4.1) a(t)= Y^[o,e)(i-jA) 

j=-oo 

is typically used. The energy of the softened pulse is 

oo 

(1-4.2) ^Q, 



]=-oo 



which tends to infinity as e ^ 0. We think of the hard pulse as an ideal non-physical pulse cor- 
responding to e = 0. In Appendix |Bl we explore the relationship between the hard and softened 
pulses. 

For a pulse of the form (I1.4.1I ). it is possible solve the Bloch equation (II. 1.11 ) explicitly. In the 
limit, as e approaches 0, the Bloch equation can be replaced by the recursion 

M_{z- (i + 1)A) = PAR^^M^{z;jA), 
10 



where Pa corresponds to a certain frequency-dependent rotation around the z-axis, and R^^^ corre- 

Re uji 

sponds to a certain frequency-independent rotation around the 



Im ujj 




-axis. The precise defini- 



tions of the operators Pa and R^. are given in Section [ 

For hard pulses, we evaluate the magnetization M_ only at the discrete time points jA for 
j e Z. It is easy to see that, at such time points, M_ is a periodic function of the frequency - it is 
a function of w = e*^^. We suppress the dependence on A and z and write M_ = M_ {w;j) as a 
function ofw^S^ and j e Z. The resulting magnetization profile is 

M{w) := lim P7^M^{w;j), 

j-»oo 

a function on . We define the discrete selective excitation transform by 

T^''^^ := M 



lim P.^M_ 



A condition on Q, that guarantees the existence of this limit is given in Section 12.51 The transform 
q-Amc jjj^pg j^^^j pulses (sequences of complex numbers) to periodic magnetization profiles (S*^ 
valued functions on S^). 

1.5. Main results: continuum theory 

The problem of inverting the selective excitation transform has been solved using inverse scat- 
tering theory. We prove in Section l2.12l 



Theorem 1.5.1. LetM. 



M 



y 



r(C) 



: M ^ 5^ C a magnetization profile, and set 



1 + M,(20 ■ 
Assume that r and are in (M). 

(a) There are infinitely many pulses uj such that Tlo = M. These pulses are parameterized by their 
bound state data (see below). 

(b) If r has the form r{S^) = e~^'^''ro(^) where tq has a meromorphic extension to the upper half 
plane with finitely many poles, and if 

lim ro(C) = 0, 

ICI^oo 

then there exists a pulse lo with rephasing time p such that Tuj = M. 



2*€p^,/or A-landB in i7+(M) (see Section^TB 



(c) Suppose that r has the form r{^) 
satisfying 

\A{0\^ + \B{Of = lforall^eR, 
and suppose that the Fourier transform of B is supported on the interval [0, 2r] C M. Also assume 
that A has an analytic extension to the entire complex plane. Then there exists a pulse uj with 
duration T and rephasing time p such that Too = M. 

The function r : M ^ C is called the reflection coefficient. The magnetization profile M is 
uniquely determined by the reflection coefficient: 

2Re r 



(1.5.1) 



miz) 



2Imr 

T+FF 
T+FF 



(-). 
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The bound state data mentioned in part (a) of Theorem 11.5.11 is defined in Section 12.31 In the 
basic and generic case, this bound state data takes the form 

D = (6)6) • • • )^m; C[,C2, ... , C^), 

where 6 ) • • • j '^m are distinct complex numbers in the upper half plane called energies, and €[, . . . ,C'„ 
are non-zero complex numbers called norming constants. The main result is that to each such r and 
D, there corresponds a unique pulse. The special case of no bound states (m = 0) corresponds to a 
pulse called the minimum energy pulse (see |^). 



1.6. Main results: discrete theory 

There is a completely analogous theory for the discrete selective excitation transform. The 
following Theorem is proved in Section |2. Ill 



Theorem 1.6.1. Let M 
and set 



M, 

my 



51 



52 c 



be a periodic magnetization profile. 



r = -. 

1 + M, 

Assume that r is in H^{S^). 

(a) There are infinitely many hard pulses Q such that T'^"'^0 = M. These hard pulses are parame- 
terized by their bound state data (see below). 

(b) If r has the form r(w) = w~PrQ(w) where tq has a meromorphic extension to the unit disk 
which vanishes at the origin, then there exists a hard pulse Q with p rephasing time steps such that 
j-dhc^ = M. 

-pB{w) 
A{w)' 



(c) Ifr has the form r{w) 



w 



y where A and B are polynomials of degree T which satisfy 

\A{w)\^ + \B{w)\'^ = Iforallwe , 

then there exists a hard pulse with duration T and p rephasing time steps such that T'^''"T2 = M. 

Again, r : S"^ ^ C is called the reflection coefficient, and the magnetization profile M can be 
obtained from the reflection coefficient using 



r 2Re r 



(1.6.1) 



M(w) 



1+ 


r 


2 


2Imr 


1+ 


r 


■2 


1- 


r 


2 


1+ 


r 


2 



The bound state data mentioned in part (a) of Theorem II. 6. II is defined in Section flAl In the 
basic and generic case, this bound state data takes the form 

D = {wi,W2, . . . ,Wm;Ci,C2, . . .,C^), 

where wi, . . . , Wm are distinct complex numbers in the unit disk called energies, and c'^ , . . . , are 
non-zero complex numbers called norming constants. The main result is that to each such r and D, 
there corresponds a unique hard pulse. The special case of no bound states (m = 0) corresponds 
to a hard pulse called the minimum energy hard pulse. In Section [2.131 we discuss the relationship 
between the continuum bound state data and the discrete bound state data. 
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1.7. Energy Formulas 

The following theorem is proved in IS. The corollary is immediate (see Section [2.11l ). 

Theorem 1.7.1. Let r be a reflection coefficient as in Theorem \1.5.1\ and let D = (^i, • • • , '^m! Ci, ^^2' • • • ' ^ 
be bound state data. Then the energy of the corresponding pulse w : M ^ C 



\u:{t)\^dt = - / log(l + |r(0|')de + 16^/mefc 

-oo ^ J —oo k—1 



Corollary 1.7.2. Let r(^) = e ^*^'^ro(i^) be a reflection coefficient as in part (b) of Theorem 
\L5.1\ and let cu be the corresponding pulse with rephasing time p. Then the energy ofoo is 

/oo A l-OO 

|cj(t)pdt = - / log(l + |r(e)|2)d^ + 16j^4/'^j4, 
-oo ^ J — oo k—l 

where ^i, . . . , are the poles ofr^ in the upper half plane, and di, . . . , dm are the corresponding 
multiplicities. 

The following theorem and its corollary are proved in Section [2. 101 

Theorem 1.7.3. Let r be a reflection coefficient as in Theorem \1.6.1\ and let D = (wi, W2, ■ ■ ■ , Wm', c'l, C2, . . . , c 
be discrete bound state data. Let ^}(t) = Yl'fL~oo ^j^it — J A) be the corresponding hard pulse. 



Then 



^ log(l+tan2^) = — / log(l + |r(e^^)|2)d0-2^1og 

)=-"00 ^ k=l 



Wk\ 



Corollary 1.7.4. Let r{w) = w^^rQ^w) be a reflection coefficient as in part (b) of Theorem 
\L6.1\ and let Q{t) = J2j=-oo '^j^if ~ 3^) corresponding pulse with p rephasing time steps. 

Then 

00 III! I 1 C'^'^ 

log(l+tan2^) = — / \og{l + \r{e'')\^)de-2Y,dk\og\wkl 

j=-oc ^ k=l 

where wi, . . . , Wm are the poles ofr^ in the unit disk, and di, . . . , dm are the corresponding multi- 
plicities. 

1.8. Algorithms 

There are three main algorithms for producing hard pulses: the SLR algorithm, the finite rephas- 
ing time algorithm, and the DIST recursion. 

1.8.1. The SLR algorithm. The SLR algorithm, discovered independently by M. Shinnar, and 
his co-workers, and R Le Roux, can only be used to design hard pulses of finite duration, as in part 
(c) of Theorem 1 1.6.11 The input to the algorithm consists of the number of rephasing time steps, p, 
and the two polynomials A and B which must satisfy 

(1.8.1) |^|2 + |5|2 = 1 

on the unit circle. The designed magnetization profile is determined using equation (11.5.11) with the 
reflection coefficient 

. ^ -pB{w) 

r{w) = w ^ , ■ 
^ ' A{w) 

The main drawback of this method is that one does not have direct control over the magnetization 
profile (more precisely, the phase of the reflection coefficient). On the other hand, one has complete 
control of the duration, which coincides with the degree of the polynomial B. The SLR technique 
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uses the fact that the magnitude of B is related to the z-component of the resulting magnetization 

' M, ' 

by the formula 



profile M 



My 



= = I^P - \B\^ = 1 - 2|Sp. 

1 + \r\^ 

The polynomial B is designed so that 1 — 2\B\'^ approximates the ideal z-magnetization, and then 
A is chosen to satisfy equation (11.8.11) . The phase of the transverse magnetization M^- + iMj, 
is not directly specified by this procedure, and ad hoc methods must be used to obtain a good 
approximation of the desired phase. This disadvantage is discussed further in Section [3.2.31 where 
we compare SLR pulses to inverse scattering pulses. 

1.8.2. The finite rephasing time algorithm. The finite rephasing time algorithm resembles the 
SLR algorithm, but it is more flexible. It produces pulses of infinite duration, and finite rephasing 
time. Similar algorithms have been presented by various authors. See for example [16, 2J. 

The finite rephasing time algorithm can be used to produce the pulse from part (b) of Theorem 
11.6.11 The input to this algorithm is the number of rephasing time steps, p, and a rational function 
ro = ^, where P and Q are polynomials (possibly of very high degree). These polynomials do not 
need to satisfy any kind of equation like (11.8.11 ). The designed magnetization profile is determined 
using equation (11.5.1! ) with the reflection coefficient 

r{w) = w '^yr^. 

Q(w) 

With this algorithm, one has direct control of the entire magnetization profile (not just the z- 
component as in the SLR algorithm), but control on the pulse duration is sacrificed. In fact, the 
designed pulse almost always has technically infinite duration. However, we will see in Chapter 
[3] that for many common applications, the effective duration of pulses designed with this method 
is within a reasonable range. The loss of direct control on the duration causes no disadvantage in 
practice. 

This algorithm is more general than the SLR algorithm, and equally efficient. In Chapter [3] 
we give a simple derivation of the algorithm, and we describe several applications in NMR pulse 
design. 

1.8.3. The DIST recursion. The discrete inverse scattering transform (DIST) algorithm, which 
we introduce in this thesis, efficiently handles the most general case of hard pulse design. It is moti- 
vated by the inverse scattering method of pulse design (see [[H). The input to the algorithm is called 
the scattering data and consists of: 

(a) An arbitrary reflection coefficient, r : —>■ C; 

(b) Arbitrary bound state data (see Section [2!4l ). 
For technical reasons, we assume that r is in H^{S^). 

The derivation of this algorithm, and the proof that the output pulse has the correct scattering 
data is the main result of this thesis. This derivation and proof can be found in Chapter [2] 
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CHAPTER 2 



Scattering theory 

2.1. Preliminaries 

2.1.1. The projection operators 11+ and n_. The Paley- Wiener theorem states that every 
function / G L^(R) can be written uniquely in the form 

/ = /+ + /-, 

where /+ G L^(M) has an analytic extension to the upper half plane, and /_ G L^(M) has an 
analytic extension to the lower half plane satisfying 

lim /±(0=0. 

We define the projection operators 11+ , n_ : L^(R) L^(M) by 

(n+/)(0 = = ^J^ me'^'dy 

(n_/)(0 = MO = ^fjiyy^'dy. 

Let and i?_(M) denote the ranges of 11+ and n_, respectively. 

Similarly, every function / G L'^{S^) can be written uniquely in the form 

/ = /o + /+ + /-, 

where /o G C is a constant, /+ G L'^{S^) has an analytic extension to the unit disk, vanishing at 
the origin, and /_ has an analytic extension to C \ D, vanishing at oo. We define the operators 

n+,n_,n+,n_ : l'^{s^) L^{s^)hy 

oo 

(n+/)H = /+H = 5^/(jK 

-1 

(n_/)(^) = Uw) = J2 /0>' 

j=-oo 

n+/ = /+ + ^/(o) 
n-/ = /- + ^/(o). 

Let i?+, and denote the ranges of 11+ , n_, n+, and n_, respectively. 

The following lemma will be needed later in the chapter. 

Lemma 2.1.1. Suppose that f is in L'^{R) and that is in L (M). Then 

(2.1.1) m+f-U+iU = ~^fiO), 

15 



and 

(2.1.2) ien_/-n_ie/ = ^/(o)- 

Proof. First note that the hypothesis implies that /(O) is well defined. We will assume suf- 
ficient regularity and decay for / so that the below integrals make sense. The general result then 
follows by continuity. 

Integrating by parts gives 

— / f'{t)e''^dt + ^ / f{t)e''^dt = -—f{0). 

Notice that the first term on the left is —H+i^f and the second term is i^n+/ . This proves (12.1.11 ). 
A similar computation can be used to prove (12.1.21) . □ 

2.1.2. The reflection coefficient. Later in the chapter we will be working with functions of the 
form r = ^ where a and b satisfy |ap + |6p = 1 on M. In this section we mention some results 
which are well known in the inverse scattering literature. 

Proposition 2.1.2. Let r G H^{R), let ^i, . . . ,^rn ^ and let di, . . . ,dm be positive inte- 
gers. Then there exist unique continuous functions a, 6 : M — > C such that 

(i) r=^on M; 

(ii) |a(OP + = 1 for all ^ G M; 

( Hi) a has an analytic extension to the upper half plane with zeros at ^i, . . . of orders di, . . . , dm,' 
(iv) limi^i^oo a(6 = 1- 

Proof. We need only consider the case where a is non-vanishing (m = 0), because adding 
zeros simply amounts to multiplying a and 6 by a common Blaschke product. Therefore, log a is an 
analytic function which, by condition (iv), must tend to at oo. We know that 

Re log a = log \a\ = —— log(l + 

Notice that log(l + |rp) is in i7^(R). Therefore, we must have 



The general solution is 



and 



loga = -n+(l + |r|2). 

m / ^ e \ dk 



ra. 



□ 



Proposition 2.1.3. Let r e H^{S^), let wi,... , Wm e and let di, . . . ,dm, be positive 
integers. Then there exist unique continuous functions a,b : ^ C such that 

(i) r = ^onS\- 

(ii) |a|2 + |6|2 = ion S\- 

(Hi) a has an analytic extension to the unit disk with zeros at wi, . . . , Wm of orders di, . . . , dm,' 
(iv) a(0) > 0. 



Proof. The proof is essentially the same as the proof of Proposition 12. 1 .21 The general solution 

™ / * \ dk 

l^\\wk\l-wlwj 



k=l 

16 



and 



ra. 



□ 



2.1.3. Banach derivatives. In this section we define tiie derivative of a curve in a Banacli 
space, and prove some results wfiicli will be needed in Section l2.12l 

Let X be a Banach space. A function n : M — > X is said to be dijferentiable at to £ 1^ if there 
exists e X such that 



lim 



ut - Ut 







t-to 



Ho 



0. 



In this case u[^^ is called the derivative of u at to- This derivative is necessarily unique. 

Let X and Y be Banach spaces, and let C{X, Y) denote the space of bounded linear maps from 
X to Y. Since C{X,Y) is a Banach space (with the operator norm), we can also speak of the 
derivative of ^ : M ^ C{X, Y). 

The following two lemmas can be proved directly by applying the above definition of derivative. 



Lemma 2.1.4. Let At -.M.^ C{X, Y) and ut : 
then so is vt = AfUt - ^ ^ Y, and 



'to 



X . If At and ut are dijferentiable at Iq, 



Lemma 2.1.5. Suppose : M ^ C{X,Y) is invertible near to £ M, and assume that both 
At and A^^ are uniformly bounded in a neighborhood oft^. If At is differentiable at to, then so is 
A~[^, and 

^)to = ~\^^to\^- 
The following propositions will be used in Section 12.1.41 

Proposition 2.1.6. Let X be a complex Hilbert space. Suppose that At ^ C{X, X) and 
f t : M — > X are differentiable at to- If has the form 

At = l + Kt, 

where Kt is a bounded, self adjoint, positive operator in a neighborhood of to, then At is invertible 
in a neighborhood of to, and the function : M ^ X given by 

Ut = A^'^vt 

is also differentiable at to- 

Proof. Since A is differentiable at to, it is certainly uniformly bounded near to. Therefore, by 
Proposition I A.O. lOl and its proof. At is invertible, and A~[^ is uniformly bounded in a neighborhood 
of to. The result now follows from Lemmas |2 . 1 .4 1 and 12.1.51 □ 

Proposition 2.1.7. Let {X, ||-||) be a Banach space of complex-valued functions on some set 
S, and suppose that there is a constant C > such that \f{p)\ < C \ \f\\for all f €z X andp €z S. If 
u : M ^ X is differentiable at to, then for each point p G S, the function 1 1— > Ut{p) is differentiable 
at to, and its derivative is given by u't^ (p). 



Proof. For simplicity, assume to = 0. We want to show that 



lim 



t 



UoiP) 



0. 



But this follows immediately from the definition of derivative, and the hypothesis that \f{p)\ < 
C7II/II forall/ G XandpG 5. □ 
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Proposition 2.1.8. Let r G i/^(K), and suppose that is in ff^(M). Then the curve t 
re*^* G (M) is dijferentiable at every to G M, and its derivative at to is given by i^re^^^° . 



Proof. Without loss of generality, we can assume that to = 0. We need to show that 



lim 



lim 



- 1 
it 



0. 



i I . One can show that there is a constant M > such that for all < t < 1 



Set MO --y 

and for all ^ G M, we have IMOl < ^ ^^d \(l)'t{C)\ < M. Given e > 0, we can choose > 
large enough so that 

||X(-oo,-c,]uK„oo)^'^|li2 > ||X(-oo -5,]uK.,oo)^5(C?')||i2 < e, 
and < t < I small enough so that 

We have 

WMrWli = \\Mr\\l2 + W^'t^r + ^td^{ir)\\l, 
(2-1-3) < ||(/>ter||i2+2||0;er||^,+2||0t95(er)||2,. 

We estimate each of the three terms in (12.1.31) : 



|X(-oo,~c.]uK.,oo)0tC?'|li2 + \\X[~ti,,^,]Mr\\l2 



< M'e' + e'Ur\\i2 



L2 



L2 



Therefore we have \\Mr\\m ^ ^(^^^ + ll^^ll^Oe^- 



□ 



2.1.4. The Marchenko equation. Later in the chapter we work with the Marchenko equation: 
(2.1.4) (1 +n+r*n_r)L = -n+r*, 

where r is a complex function on A = M or 5^. The Marchenko equation comes from a system of 
equations: 

K = 1 + U+r*L* 
L = -U+r*K*. 

In the literature A is M, and the system is typically written in the Fourier domain: 



HK-l)ix) 
L{x) 



fi-x-yTLiyTdy 

oo 

fi-x-yykiyydy, 



where / = f. In this section, we discuss hypotheses on r which guarantee that there is a unique 
solution L G H]_{A) to equation ( 12.1.41 ). 

Lemma 2.1.9. //r G i^^A), then A := n+r*n_rn+ is a bounded, positive, self adjoint 
operator from (A) to itself. 
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Proof. By Fact IA.0.41 multiplication by r or r*is a bounded operator from H^{A) to itself. 
It is clear that n_ and 11+ are also bounded operators on H^{A). Therefore A is itself such an 
operator. We can use Fact lA.O.Tj to show that A is positive and self adjoint once we establish that 
(n+r*n_, n_rn+) forms an adjoint pair. This can easily be shown by proving that (n_, 11+) and 
(r, r*) are each adjoint pairs. □ 

Proposition 2.1.10. If r £ H^W, then there is a unique solution L £ ^^+(A) to the 
Marchenko equation H2.1.4\l . The norm of this solution satisfies the estimate 



\L\\m < 



Proof. The Marchenko equation can be written 

(1 + A)L = -n+r* 

for A defined in Lemma 12.1.91 The Lemma tells us that ^ is a bounded, positive, self adjoint 
operator from (A) to itself. The desired result then follows immediately from Proposition lA.O. lOl 



□ 



Proposition 2.1.11. Letro G H^{M), and let A : M C{H^{M.), H^{M.)) be given by 

At = n+rSe-2*?*n_roe2^«*. 
T/'^ro is in i^^(M), then the solution L : M ^ if^(R) to the Marchenko equation 

{l + At)Lt 

is dijferentiable (in the sense of Section \2.1.3\l at every to G 



Proof. By Proposition 12.1.61 we just need to show that A is differentiable at to- It is suffi- 
cient to show that multiplication by tqc^*^* is differentiable at to- This can easily be shown using 
Proposition |2X11 □ 

2.2. The Zakharov-Shabat system 

The selective excitation transform has been shown to be equivalent to the scattering transform 
for the Zakharov-Shabat (ZS) system of equations. In this section we give the details of this rela- 
tionship. We mainly follow the notation of |5|. 

The Bloch equation (II. 1.11 ) can be written as 

z — Im Lo{t) 

dtM{z;t) = -z Rea;(t) M{z;t). 

Imcj(t) -Reu;(t) 

We can think of the 3x3 matrix in this equation as an element of the Lie algebra of 5O3M mapping 
M G 5^ to an element of the tangent space TmS'^. If we lift to the universal cover, SU2C, this 
equation becomes 



(2.2.1) 
where 

(2.2.2) 

and 
(2.2.3) 



dt 







-ii q{t) ' 








-q%t) ii 








2Re il)lij2 




M{z;t) 




2Im 7/^5-^2 













Qit) 



-y{t). 



The function q is called the potential for the ZS-system. 
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to the differential equation 



-ii q{t) 
-q*{t) 



2.3. Continuum theory 

In this section we outline the scattering theory for the ZS-system. Many of the formulas can be 
found in (5). 

Let g : M ^ C be an integrable potential for the ZS-system. Then there exist solutions i/ji± = 

(2.3.1) 
satisfying 
(2.3.2) 
The matrix 
(2.3.3) 

is independent of t, and is called the scattering matrix. Let us define A± and B± by 



lim e*«Vi±(e;t) 

t— >±oo 



for all i E 



(2.3.4) 
so we have 
(2.3.5) 
and 
(2.3.6) 



lim 

i— >itoo 







" 1 " 










for all ^ E 



AO 



AXAO BiAO 
-B+AO A+AO 



A^AO -BiAO 
B^AO A*iO 



for all t and ^. One can show that the resulting magnetization profile, M = Tlo, from Section [TTT] 
is given by equation ( 11.5.11 ) for the reflection coefficient 



AO 



m 

AO' 



The following is an outline of the main elements of the scattering theory (see [^."Sl). The Marchenko 
equations below do not appear in their typical forms. The derivations of these Marchenko equations 
are given in Sections |2!6] and f2J\ 

• The functions a and b satisfy 

i«(e)i' + iK6i' = i foreeffi 



(2.3.7) 



(2.3.8) 



(2.3.9) 



" AO ■ 




" 1 " 


. biO . 








and 

lim 

For every t, ^ E M, we have 

\A±AO\' + \B±AO\' = i- 

For each t, the functions ^4^ ^, B*^ ^, A^^t, and B^ t have analytic extensions to the upper 
half plane. 
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The function a = A^^^A^^t + has an analytic extension to the upper half plane. 

We assume that a has finitely many zeros {,^1, . . . , ^„i} in the upper half plane, which are 
all simple. For each zero ^fc, of o, there is a constant C{. such that 



(2.3.10) 



(2.3.11) 



B-Mk) 



for all t e 



Set 



and 



(2.3.12) 



Ck 



a'(6 



a' (6) Ck[a'{ik)?' 

The data (a, 6; ^1, ... , C[, . . . , C!^) is called the scattering data for the potential q. 
The data (r; ^1, . . . , C*!, • • • , Cm) is called the reduced scattering data for the poten- 
tial q. The functions a and b can be determined from the reduced scattering data by the 
formulas 



n 

k=l 

ra. 



exp(-n+(l + |r|2)) 



• The function i?^ ^ can be determined from the reduced scattering data. It is the unique 
solution (see Proposition 12. 1 . lOl ) to the Marchenko equation: 



(2.3.13) 



where 



(2.3.14) 



(i + n+r;n_rOi?;,i = -n+r,* 



rt(^) = n_re 



fc=i 



(2.3.15) 



• The function ^ can be determined from the left reduced scattering data 

S = (s; ^1, . . . , £,m\ Cl, . . . , Cm), 

where 

6*(0 



«(0 



It is the unique solution (see Proposition 12.1.101) to the left Marchenko equation: 
(2.3.16) (1 + U+s:U_st)B_^t = -n+s,* 

where 



(2.3.17) 



k=i 



(2.3.18) 



The potential q can be recovered using 

1 



or 



(2.3.19) 



q{t) = -HBI,){0' 



q*{t) = -HB-,t)iO^). 

TT 



The following is a restatement of part (a) of Theorem [T37T] in terms of the ZS-system framework. 
The proof is given in Section 1X121 
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Theorem 2.3.1. Let S = (a, 6; ^i, . . . , . . . , C^) be arbitrary scattering data, as 

above, such that r = ^ and are both in H^{M). Then there is a well defined potential q for 
the ZS-system such that S is the corresponding scattering data. This potential can be found ei- 
ther by using equations M.S. 13^ , M.S. 14^ and M.S. 18^ . or by using equations M.S. 16^ , M. S.17i and 



Remark 2.3.2. Typically the right Marchenko equation ( 12.3.131 ) is used for the positive values 
of t, and the left Marchenko equation (12.3.161) is used for the negative values of t. 



2.4. Discrete Theory 

In this section we describe an analogous scattering theory for hard pulses. 
Consider a potential of the form 



Q{t) = ^ fikS{t-j6), 



j=-oo 



such that 



j=-oo 



We will call such a function a discrete potential. For these potentials, the differential equation 
(12.3.11 ) is replaced by a recursion: 



(2.4.1) Vi±(e;(i + i)A) 



e-'^^ 
e^^« 



cos I flj I 

■ I 

■ sm 



cos I 



V'i±(C;jA). 



|^oiii|//j| COS I I 

For each integer j, (,^; j A) and ^2- J A) are periodic functions ofw2= e'^^. Let us set 



(2.4.2) 

Then the recursion is 
(2.4.3) ^'±,i+i(w^) 
and the scattering matrix is 
(2.4.4) 
Let us define 
(2.4.5) 

so we have 
(2.4.6) 
and 



V'i±(C;jA). 



_ 1 

W 2 



Wi 



cos Ifljl 



'a -b* ' 












b a* 













'2-,j 



A±j{w)w 2 



lim 







' 1 " 


B±,,{w) 








for all w e S"" 



-b*w-^ 
a* 



A* R* 



-B* 
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for all j G Z. One can show that the resulting magnetization profile from Section [Ol is given by 
equation (11.6.11) for the reflection coefficient 



r{w) 



b{w) 
a{w) 



The following is an outline of the main elements of the discrete scattering theory. Some of these 
statements are proved in Section 1231 The discrete Marchenko equations are derived in Section |Z81 

• The functions a and b are in L^{S^) and satisfy 

laHp + = 1 for weS^ 



(2.4.8) 



and 



(2.4.9) a(0) > 0. 

• The functions A± and B± are in and satisfy 

(2.4.10) \A±j{w)f + |S±j(u;)|2 = 1 for weS^ 
and 

(2.4.11) i±j(0)>0. 

• For each j, the functions A*^ ^, B"^ ^, A^^t, and w^^B^^t have analytic extensions to the 
unit w-disk D. 

• The function a = A^jA^j + B*^jB^j has an analytic extension to the unit disk. We 
assume that a has finitely many zeros {wi, . . . , Wm} in the unit disk, which are all simple. 
For each zero Wk of a, there is a constant c'^ such that 



(2.4.12) 



^-j(wfc) 
B-,j{wk) 



4 



-B*+j{wk)wl 
A*{wk)wl 



for all j eZ. 



Set 



(2.4.13) 



and 



(2.4.14) 



Cfc 



a'{wk) 



-w. 



a'{wk) 



Ck [a'{wk)f 

• The data {a,b;wi, . . . , Wm', c'l, • • • , c^) is called the discrete scattering data for the po- 
tential q. 

• The data {r;wi, . . . , Wm', ci, . . . , Cm) is called the reduced discrete scattering data for the 
potential q. The functions a and b can be determined from the reduced scattering data by 
the formulas 



(2.4.15) 



n 

k=l 

ra. 



Wk - w 
\wk\ 1 - wlw 



•exp(-n+(l + |rp)) 



The function 



wB 



can be determined from the reduced scattering data. It is the unique 



solution (see Proposition 12. 1 . lOl) to the Marchenko equation: 



(2.4.16) 



(i + n+r*n_r,; 
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where 



(2.4.17) 



rj = U^rw^ ^ — 

k=l 



W -Wk' 



The function 
where 



A-.t(O) 



can be determined from the left reduced scattering data {siwi, . . . ,Wn;ci, 



(2.4.18) 



a 



(2.4.19) 



It is the unique solution (see Proposition 12.1.101) to the left Marchenko equation: 



(i + n+s*n_sj 



^-,i(0) 



where 



(2.4.20) 



k=l 



W-Wk 



The potential q{t) = ^j5{t — j A) can be recovered using 



(2.4.21) 



l7,l 



arctan |7j| 



for 



(2.4.22) 



wB*.- 
lj=T{ . +-^ (1) 



or 



(2.4.23) 



7.: = ^i- 



5_ 



■)(!)• 



^-j+i(0)- 

The functions and -B+j can also be approximately computed recursively. We start 
by setting 



(2.4.24) 





-,Af(^«) 




" 1 " 












for M >> 0, and then use the recursion 



(2.4.25) 



B. 



for 



(2.4.26) 



7,- 



:i + i7,i 



l7. 



-7i 



(2.4.27) 



Similarly, the functions A-^j and can be approximately computed recursively. We 
start by setting 

B^^^m{w) 



for M >> 0, and then use the recursion 



(2.4.28) 



A. 
B_ 



(l + l7,P)-^ 
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' 1 " 







1 


_ -1 





7i 



A. 
B. 



for 



(2.4.29) 



7j 



The recursions (12.4.251 ) and (12.4.281 ) along with equations (12.4.261 ) and (12.4.291 ) form the discrete 
inverse scattering transform (DIST) algorithm. These equations are derived in Section [Z91 

The following is a restatement of part (a) of Theorem 11.6. II in terms of the ZS -system frame- 
work. The proof is given in Section 12.111 

Theorem 2.4.1. Let S = {a,b;wi, . . . , Wm', ci, . . . , Cm) be arbitrary discrete scattering data, 
as above, such that r = ^ is in H^{S^). Then there is a well defined discrete potential q{t) = 
Yl'^-oo fJ'j^it ~ 3^) for the ZS-system such that S is the corresponding discrete scattering data. 
This potential can be found either by using equations l \2.4.16i , A2.4.17i and H2.4.22i . or by using 
equations <\2.4.19l <\2.4.20\l and A2.4.231 



2,5. Forward discrete scattering 

In this section we prove the analyticity properties of A± and B± from Section [2!4l 
Proposition 2.5.1. Let ^ -."L ^ Cbe a sequence of complex numbers such that 



< CXD. 



Then there are unique solutions A± and B± to equations l \2.4.25i . l \2.4.28i and H2.4.6i . Further- 
more, for each integer j, the functions A^j, w~^B-j, A"^ ■ and B*^_ ■ are all in H^(S^). 



Proof. Let Aq and Bq be the solutions to the recursion 
(2.5.1) 



normalized by 



Oj + l 



Oj + l 



(l + l7,f 



1 



A 



" ^0,0 ' 




' 1 " 










Notice that equation (12.5.11) is equivalent to (12.4.281 ). Clearly we have |^ojP + I-SqjP = 1 on 5^ 
for all integers j. Therefore, we can estimate 



\Aoj+i-Ao,j\ <hj\ + 



and 



w 



which implies that the sequences Aq^, Aq^i,Aq^2, ■ ■ ■ and Bq q, w^^Bq^i, w^^Bq^2, ■ ■ ■ converge in 
L°°( 5^) to some functions ap and 5o> respectively, in (5^) C L'^{S^). One can show inductively 
that Aqj and Bq are in ^+(5'^) for all j > 0. Thus, uq and also must be in H^{S^). By 



multiplying the matrix recursion on the right by 
equal to 



h* 



'0 



, we see that A+ j and i?+ j must be 



A+j = UqAqj + bow^BQj 
B+j = UqW^^Boj 
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which impHes that ^^J, = ^0 and B*^_ q = —b^ are in ^+(5'^), as desired. By similar reasoning, 
A\j and B^j must be in H+{S^), for all j. □ 

2.6. Derivation of the right Marchenko equation 

In the next two sections we derive the right and left Marchenko equations for the ZS-system. 
We assume that q is an integrable potential with scattering data 

5 = (a, 6; 1^1, ... , .^mi Ci, . . . , C^). 

Let A± and B± be as in Section 1231 Recall that A*^_^, B*^^, A^^t, and B^ t all have analytic 
extensions to the upper half ^-plane. Rearranging equation (12.3.61) gives 



(2.6.1) 



B+AO AXAO 



AO 



or 



-B^^t = B+^t + re^'^'Al 



We apply n_, conjugate the second equation, and rearrange: 
(2.6.2) A+t = 1 + n_re^'«*5l^ + n_-A„ 



(2.6.3) 



B 



-n+r*e 



Using the properties of , . . . , we get 



N 



a 



j^Ce^^^^^BXAik) 1 



k=l 

N 



Y^Cke^^^^'BXA^k)-j^ 

^-Q+,tBl 



and 



a 



n' 



^ C'.e'^^-'AXACk) 1 



a'(6) e - 6 



k=l 
N 



k=l 

^-Q+,tA\^t 
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where 



Q+AO = E 



k=l 



Therefore (I2.6.2I ) and (12.6.31) become 

(2.6.4) A+^t = l + n_(re2^«*-Q+,t)5;^ 

(2.6.5) = -n+(re2^«*-Q+,t)*^+,t 
or 

(2.6.6) A+^t = l + U_rtBl, 

(2.6.7) 5;^, = -n+r;A+,t 

where 

n(0 = n_re2*«*-Q+,i(e) 



k=l 



Equations (12.6.61) and (12.6.7b can be combined into the single Marchenko equation 
(2.6.8) (1 + n+r*^U_rt)Blt = -U+r^ 

To prove equation (12.3.181) we rewrite the ZS-system in terms of A± and B±: 



qit) 
-q*{t) 2ii 

Taking the f-derivative of both sides of equation (12.6.71 ) gives 



d 






Ik 







q{t)A* - 2iiB\ 



d 



-n+(-ri)M+,t-<?(t)n+r;i?+,t 
Ii+2iirlA+^t-q{t)AXt+q{t). 



Here we used the fact that 



dt 



We then use equation (12.1.11 ) from Lemma l2. l.ll to solve for the potential: 

q{t) = 2ien+r,M+,t - n+2ier?A+,i 



1 



1 



TT 



•^(^;,t)(o+). 



2.7. Derivation of the left Marchenko equation 

A similar method can be used to derive the left Marchenko equation. Instead of equation ( 12.6.11 ). 



we use 



or 



a{0 



-A 



a 
1 



-B. 



b* 



B*_ 



-2i. 
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Again, we apply n_ and conjugate the second equation: 

(2.7.1) A*_t = l + n„se~^^«*B_j + n_i^* , 

a ' 

(2.7.2) B^t = -U+s*e^'^^A*_f + U+—B+t. 

a* 

This time, we have 

" AU(t) 1 



and 



-n_g_,tS_,i 



where 



fc=l 

Therefore (12.7.11 ) and (12.7.21) become 



fc=i ^ ^'^ 



= l + n_(se-2'?*-Q_,j)S_,j 
= -n+(se-2*€*-Q_,i)Ml,, 

or 

(2.7.3) A*_^t = l + U^stB^^t 

(2.7.4) = -n+sJ'Al^t 
where 

stiO = n_^e-2*«*-Q_,t(e) 

6* ^ C.e-^^?^* 



Equations (12.7.31) and (12.7.41) can be combined into the single Marchenko equation: 

(i + n+Si*n_st)s_,j = -n+s*. 

The proof of equation (12.3.191 ) is identical to the above proof of equation (12.3.181 ): 
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-q*{t)A.^t + 2iiB.,t = -n+(-dt)M*__, - q*{t)Xl+s*B*_^^ 

= -n+2ieSiM*_,-g*(t)A_,t + g*(t) 

q*{t) = -2iCU+s:A*_^,+U+2iCs:A*_^, 

= i^(.*A*_,)(0) 
vr 

= -i^(B_,,)(0+). 
vr 

2.8. Derivation of tlie discrete Marchenko equations 

In this section we derive the right Marchenko equations for hard pulses. We omit the derivation 
of the left equation, but the reader should be able to reproduce it using the techniques from this 
section and the previous two sections. 

Assume that q has the form 

oo 
j=-oo 

where 

oo 

j=-oo 

and let S = {a,b;wi, . . . , Wm'-, c'l, • • • , c^) be the corresponding discrete scattering data. Let A± 
and B± be as in Section lZ4l Recall th at B'^ A^j, and w^^B^j all have analytic extensions 
to unit if-disk. Rearranging equation (12.4.71 ) gives 

bw^ a* 
or 

= A+ j - rw^-\wBl .) 
-w'^B_i = w-^B+i +rw^-^Ali. 

We apply n_ and conjugate the second equation: 

(2.8.1) A+j = A+j{0) +U^rw^-\wBlj) +U^^A^j 

(2.8.2) wBlj = -U+{rw^~^yA+j +U+^wB*_j. 
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A. 
B. 



-B* 
A* 



A, 
B 



Using the properties of c'^, . . . , c'^, we get 



N 



A-j{wk) 1 



k=l 

N 



a'{wk) w-Wk 



k=l 



W -Wk 



and 



a 



^ a'{wk) w -Wk 
j^d,nP-'AX^{wk) 1 

k=l 
N 



a'{wk) w-Wk 
y2ckwi~^Alj{wk) 



k=l 



W-Wk 



where 



m j — l 



k=l 



W-Wk 



So, equations (I2.8.1I ) and (12.8.21) become 



(2.8.3) 
(2.8.4) 
or 

(2.8.5) 
(2.8.6) 

where 



A, 

wB: 



i+,,(0) + n_(™^-i - Q+^,){wBl^^) 

-n+(™^-i-Q+,,)M+j 



A 



= A+^,{G) + Ii^T,{wB\^^) 



wBXj 



= -U+r*A+j 
U_rw^^^ - Q+j 

N j-1 

TT 7-1 ^kW% 



k=l 



W-Wk 



Equations (12. 8. 5b and ( 12. 8. 6b can be combined into the single Marchenko equation: 
(2.8.7) (l + n+r*n_rj) 



wBl ■ 



To prove equations (12.4.221 ) and (12.4.231 ). we write the recursion (12.4.11 ) in terms of A± and B± : 



(2.8.8) 



B±,j+i 



(i + i7,fr^ 
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- 1 7/ 




' A±,j ' 


—JjW w 







(2.8.9) 
for 



(l + l7,l 



These recursions immediately imply 



and 



as desired. 



1 



w ^ 



7, = ^tan|;x,|. 



B 



-77 



(1) 



(0) 



2.9. Derivation of tlie DIST recursion 

To derive equation (12.4.221) we simply plug equation (12.8.61) into the recursion (12.8.' 



wBl 



Examining the coefficient of u;^, we have 

7, (^(^+,,+i)(0) - :F{w-'r*{w-'B+j+^)m) = V*^+,,-+i)(0). 

This immediately gives the desired result. 

A similar computation can be used to obtain equation ( 12.4.231 ): 



—jjwA^j + wB^j 
-j*{wA^,-U+s*^,Bl^^) 

-7* {HA-,m-Hw''^UiB*-,jm) 



U+s*^,A*_^^-j*U+s*^,Bl^^ 
wB^ o 



^* ,)(o). 



2.10. The discrete energy formula 

In this section, we prove Theorem 11.7.31 and Corollary 1 1.7. 41 Let 

S = {a,b;wi, . . . ,Wm;c[, . . . ,c'^) 

be discrete scattering data, and let A± and B± be as in Section [2!4l By equation (12.4.71) . we know 
that a is given by 

a = lim A_ 



Therefore, since limj^_oo A^ 



1, the recursion (12.8.81) tells us that 



a(0)=a(0)= n (1 + 17.1 



Let us write 



a{w) = ao{w) ■ 



W — Wk 

L \wk\ 1 - wlw' 
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where oq is analytic in the unit disk, and wi, . . . , Wm are the zeros of a in the unit disk. Since 
log |ao| is harmonic in the unit disk, and since a(0) is positive, we have 

f27r 



^ / \og\ao{e'')\de = log|ao(0)| 
27r Jo 

m 

= loga(O) - E^ogl^fcl 



k=l 

oo 



E log{l + \-fj\^)-f2^og\wk\. 

j=—oo k=l 



Therefore, 



oo 27r "I 

(2.10.1) E + = -- / ^og\ao{e'')\de-2^log\wk\ 

A- ^ JO f,^l 

1 



For ^t; on the unit circle, we have |ao(w)| = |a(t(;)| = (l + |r(z/;)p) 2,wherer = Ms the reflection 
coefficient. Therefore ( 12.10.11 ) becomes 



E log(l + |7,f ) = ^ / log(l + |r(e^^)|2)d0-2Elog|^fc| 

j=-oo '—1 

By equation (12.2.31 ). this is 



k=X 



oo 111 /'27r 

E log(l+tan2^) = — / log(l + |r(e*^)|2)(i0-2Elog|«^A:|. 

This proves Theorem |1.7.3l Corollary 11.7.41 follows immediately in light of the proof of part (b) of 
Theorem ll.6.11 Notice that if the time step. A, is small, then the magnitude of ujj is also small, and 
so we have 

oo . „27r ™ 

E l^il'^^/ log(l + |r(e^^)|2)d0-8Elogl^A:|. 
i=-oo ^ 1^^^ 

2.11. Proof of the main result: discrete case 

In this section, we prove Theorem ll.6.11 

Proof of part (a): We start with the reduced discrete scattering data 

S = {r;wi, . . . ,Wm;ci, ...,Cm) 

where r is in H^{S^). This corresponds to unique discrete scattering data 

S = {a,b;wi, . . .,Wm;c[, . . . ,c'^). 

Proposition 12 . 1 1 . 6l belo w tells us that there is a unique hard pulse with scattering data 5*. 

Proof of part (b): This follows from part (a), in the special case where, in the notation of 
Section |Z81 Q+,p = vq, because in this case we have rj = for all j > p. This implies that the 
hard pulse vanishes for time steps j > p, as desired. See Appendix ICl for the case of non-simple 
poles. 

Proof of part (c): This follows from part (b). Simply note that the scattering data in this case 
has the form 

S = {A,W'fB;Wi, . . .,Wm]Ci, . . .,Cm). 

One can check that the scattering data for the time reversed pulse pj = —p*_j is 

S = (A, -wPB*;wi, . . .,Wm;ci, . . .,Cm)- 
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Notice that w^~''w'' B* is analytic in the upper half plane, which implies that the time reversed 
pulse ends at time step T — p, as desired. We omit some details about the bound state data which 
need to be worked out. 



Lemma 2. 11.1. Let S = {a,b;wi, . . . ,Wm',ci, 



,Cm) be discrete scattering data such that 



^ is in H^{S^). For each j there exist unique functions A^j, w ^-B_j, ^+ and B 



in 



H^{S ) such that 
(2.11.1) 

(2.11.2) 



B. 



-B*_ 
A* \ 



A, 
B. 



A* 



2tt 



a 



a* 



(2.11.3) 

and 

(2.11.4) 

(2.11.5) 

where 



i±,i(0) > 0, 



a 



n. 



-w 



'B^ 



U^Q+jA 



Q+,jiw) 



E 

k=l 



CkWi 



W -Wk 



Proof. We first prove uniqueness. From Section 
solution to equation (12.8.71 ) where 

rj = Il-rw-^~^ — Q 
Using equation (12.8.51) we can determine 



wBl . 

we see that — 



(0) 



must be the unique 



A 



1 + n_r,- 



wB 



The functions A^j and B^j are then uniquely defined by (12.1 1.21 ) and ( 12.11.31 ). Finally, A_j 
and B^j can be computed from A^j and B^j using the matrix equation (12.11.11 ). This proves 
uniqueness. 



To prove existence, we just need to show that A^j, B^j, A*^ and B*^ ■ as defined above 
are all in H\{S'^), and satisfy equations (12.11.41 ) and (12.1 1.51 ). By construction, we know that A*^ ^ 

and B*^ j are in H^iS^)- Also, by construction, we know that equations (12.8.11 ). (12.8.21) . (12.8.31) . and 
(12.8.41) hold. Comparing these equations immediately gives (12.11.41) and (12.1 1.51 ). So we just need to 
show that A^j and w^^B^j are analytic in the unit disk. But this actually follows from equations 
(12.11.41) and (12.11.51 ). Indeed, equation (12.11.41 ) imphes that ^A-j + Q+jwB'^ j is analytic in the 
unit disk. Since ^ and Q+j have the same poles and multiplicities, and since wB*^ ■ is analytic 
in the disk, it follows that A_ j must also be analytic in the unit disk. Similarly, equation (12.1 1.51 ) 
proves the analyticity of w^^B^j. □ 

Lemma 2.11.2. The functions A±j and B±j given in Lemma \2.11.1\ satisfy the recursion 



A 



B± 



(l + l7jl 



- 1 7, ■ 






— l*jW w 
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where 



7,- 



^+,i(0) 

Proof. By uniqueness of the functions we just need to show that A-t j+i and given by 

A 



and 



7i - 



- 1 7/ 










. B±,j _ 



have the desired properties at step j + 1. By multiplying equation ( |2.11.1| | on the left by (1 + 



1 



7j 



it is easy to see that the matrix equation holds. It is clear that A_ j+i, 



j+i, A*^ jj^i, and B'X ^j^^ are in H^{S^), because 7^ was chosen precisely so that this 
condition would hold. Also it is clear that properties (12.11.21 ) and (12.11.31 ) are satisfied. So we just 
need to verify equations (12.11.41 ) and (12.11.51) at step j + 1: 



(1 + |7,f )^n„(ii_,,+i + Q+,,+i7x;^;,^.+i) 



,1 



(1 + h,\^)mdw-'B^,+i - 





n_(-7;iA_,, + is_,, 

^ J n '-^ '-^ 



-n-7;Q+,j+i5;,i 
,1 



0. 



Here we used the fact that 



where g is analytic in the unit disk. 

Lemma 2.11.3. The functions A^j and B^j given in Lemma \2.11.1\ satisfy 



□ 



1 



Proof. Since 



wB* 



A 

B4 



is the solution to 



(1 + U+r*U^rj 



as j +00. 



we can use Proposition 12 . 1 . 1 Ol to estimate its ff^-norm by 

wB* 



'+,3 
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For i > this is 



B 



Clearly the expression on the right tends to zero, which implies that 

B* 



lim 

j-*+oo 



0. 



By equation (I2.8.5I) and Fact lA.0.4l we have (for j > 0) 



A 



wB 



wB* 



< 2 



B: ~ 



which implies that 



as desired. 



lim 



A 



0, 



□ 



Lemma 2.11.4. Let A± and B± be the functions given in Lemma {2.11. h and let be as in 
Lemma \2.11.2\ Then for each j we have 



(2.11.6) 



A, 



+ 



B 



^+,i(0) 



n(i+i7fci')-'o« s\ 

k=j 



Proof. The recursion in Lemma l2.11.2l implies that 

i+,,+i(0) = (l + |7,f)-^i+,,(0) 

and 

+ = |A+,,f + |i?+,,f on S\ 

Together, these imply that 



A 



+ 



B, 



^+j+i(0) 



(l + l7,f 



+ 



B 



on 



for all j. Therefore, formula (12.11.61 ) follows from Lemma l2.1 1.31 

Lemma 2.11.5. The functions A^j and Bj^j given in Lemma \2.11.1\ satisfy 



□ 



(2.11.7) 

and 

(2.11.8) 



A+j - 1 



B 








as j +0O 



|yl+jf + |S+jf = Ion S^. 
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Proof. Lemma [2.1 1.41 implies that + is constant on S^. The recursion in 

Lemma (I2.n.2l ) implies that this function is independent of j. So equation (I2.11.8I ) follows from 
property (I2.11.2I ). 

Lemma |2.1 1.41 also implies that limj^oo ^+,j(0) = 1- Therefore equation (12. 11 .71 ) follows 
from Lemma 12. 11.^ □ 



Each of the above lemmas has an analogue for the functions A^j and B^j, which we omit. 
The above lemmas together with their analogues, give the following 



Proposition 2.1 1.6. Let S = {a,b;wi, . . . , Wm', ci, . . . ,Cm) be discrete scattering data such 

^^B^j, A*^ j, and 



that r = ^ is in H^{S^). For each j G Z there exist unique functions A^^j, w 
B"^ j in H]^{S^) such that 

(2.11.9) 



a 



a* 



A* 



-B~ 



A+A 



A. 
B. 



-B*_ 
A* ] 



and 

(2.11.10) 

(2.11.11) 
where 

(2.11.12) 

Furthermore, if we set 
or 

then these functions satisfy 



n_-A_ 

a 

li^-w-^B^ 



Q+Aw) 



CkW^ 



E 

k=l 



W -Wk 



7i 



-7,- 



B±,j+i 



and 



(i + i7,pr 



■ 1 7i ■ 




' A±,j ' 


^jW w 




. B±,j _ 



as J 



iboo. 



2.12. Proof of the main result: continuum case 

In this section, we prove Theorem 1 1.5 .11 

Proof of part (a): We start with the reduced scattering data 

S = (r ; ^1 , . . . , ; C*! ) • • • 1 Cm ) 

where r and are in iJ^(IR). This corresponds to unique scattering data 

S = (a, 6; ^1, ... , i^m,; Ci, . . . , C^). 

Proposition 12. 12.5l below tells us that there is a unique pulse with scattering data S. 

Proof of part (b): This follows from part (a), in the special case where, in the notation of 
Section [Z6l Q+,p = ?'o, because in this case we have rt = for all t > p. This implies that the 
pulse vanishes for times t > p, as desired. See Appendix O for the case of non-simple poles. 
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Proof of part (c): This follows from part (b). Simply note that the scattering data in this case 
has the form S = {A, e~'^^f^B; ^i, . . . , Ci, . . . , Cm)- One can check that the scattering data for 
the time reversed pulse ^{t) = —q{—t)* is 

S = {A, —e '^^B ; ^1, . . . , Ci, . . . , Cm)- 

Notice that e'^^^'^'^^^w'^^^^B* is analytic in the upper half plane, which implies that the time reversed 
pulse ends at time T — p, as desired. We omit some details about the bound state data which need 
to be worked out. 

Lemma 2.12.1. Let S = {a,b;^i, . . . ,£,m',Ci, . . . , Cm) be scattering data such that r = ^ 
and are in H^{W). For each t G M there exist unique functions A^^t — h ^+ t> '^^^ t 

in H\_ (R) such that 



(2.12.1) 

and 

(2.12.2) 



A. 
B. 



-B*_ 
A* ' 



a 



U_-A. 



-U_Q+,tB* 



a* 



(2.12.3) 
where 



n_-s_ 

a 



Q+A^) 



k=i 



Cke 



we see that B\_ ^ must be the unique 



Q 



Proof. We first prove uniqueness. From Section | 
solution to equation (12.8.71 ) where 

Using equation (12.6.61) we can determine 

A+^t = 1 + n_rtB;_,. 

Finally, ^_ j and B^ t can be computed from t and i?+ t using the matrix equation ( 12.12.11 ). 
This proves uniqueness. 

To prove existence, we just need to show that A_ t — 1, A*^ t — 1, and B*^ j as defined 
in the previous paragraph are all in i7^(M), and satisfy equations (12.12.11) . (12.12.21) and (12.12.31 ). 
By construction, we know that A\ ^ — 1 and B*j_ ^ are in //^(R). Also, by construction, we know 
that equations (12.6.21) . (12.6.31) . (12.6.41 ). and (12.6.51 ) hold. Comparing these equations immediately 
gives (12.12.21) and (12.12.31) . So we just need to show that A^^t and S_ ^ are analytic in the upper 
half plane. But this actually follows from equations ( 12.12.21 ) and ( 12.12.31 ). Indeed, equation ( 12. 12.2b 
implies that \A-^t + Q+,tB'^ ^ is analytic in the upper half plane. Since ^ and Q+^t have the same 
poles and multiplicities, and since B'^^ is analytic in the upper half plane, it follows that A_^t 
must also be analytic in the upper half plane. Similarly, equation (12.12.31 ) proves the analyticity of 
B-,t. □ 

Lemma 2.12.2. Let A±^t cind B±^t be the functions given in Lemma \2.12.1\ Then for each 
^ G M, ^±,t('C) -S±,t(0 '^^^ dijferentiable with respect to t, and we have 



(2.12.4) 
where 



dt 







" 7t ' 




" A+AO ' 


_ B+AO . 




. -l*t _ 




. B+AO . 



7t = -.F(s;,)(o+). 

TT 
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Proof. Let us first show that for each ^ e R, t(C) and i?+ t(^) are differentiable with 
respect to t. By Proposition 12.1 .71 it is enough to show that yl+ ^ — 1 and i?+ t are differentiable as 
curves in the Banach space (M). Let us prove this for Bj^ t using the fact that it is a solution to the 
Marchenko equation ( I2.6.8I ). By Proposition |2. 1.6| and Lemma [2.1.4[ it is sufficient to show that (1 + 
7r_|_r^7r_r() and vr+r^ are differentiable as curves in C{H^ {M) , (M)) and i/^(M), respectively. 
But this follows easily from Proposition 12.1.81 and Fact lA.0.4l Finally, the differentiability of t 
is apparent using equation (12.6.61 ). 

So we know that A+^j and Bj^^t have t-derivatives t and B+.t in Differentiating 
equations (12.6.61 ) and (12.6.71) with respect to t gives 



(2.12.5) 
(2.12.6) 



A+^t = Ti-2iirtBX^t + Ii^rtBX^t 

B+^t = -n_2ien^;,t-n_ni;,,. 



This system has a unique solution for j and B+^t since it can be combined into the single 
Marchenko type equation 

(1 + n_nn+r;)5+,t = -n_2ierM;,t + n_nn+2ier;s+,t 

(see Section [2T41 ). Plugging ITTZ4] into [TTTS] and [TTTe] gives 

^tB+,t = Ii-2iirtB\^^-^t'^_rtAX^-Ii_2iinBl^^ 
-jtA+,t + 2i^B+,t = -n_2i^nAl,-j;U^rtBX,. 
Using (12.6.61 ) and (I2.6.7I ). these equations become 

= 

2i^B+^t = -n_2iert^;,t + 7;- 
The second of these equations is satisfied if we set 



-It 



-2iiB+^t - li-2iirtA\ ^ 
2i£Ii_rtA\^t - Ii_2iirtAXt 

^r{nA%^,m 



1 



vr 



(see Lemme lZ 1.11 ). 

Lemma 2.12.3. The functions t and -B+,t given in Lemma \2.12.I\ satisfy 

i?^(M) 



□ 



A+,t - 1 








as t 



+00. 



Proof. Since B'^ ^ is the solution to 

(i + n+r;n_ri)i?;,t 

we can use l2.1.1(3l to estimate its i^^-norm by 



For t > this is 



\K,t\\m ^ 



Clearly the expression on the right tends to zero, which implies that 

0. 



lim 1 1 SI 

t— >+oo 
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By equation (IA.0.41 ) we have (for t > 0) 



\^+,t-l\\m 



which implies that 



^ 2 llrollji^i 



lim ||^+ i - 1|| ffi =0, 

t~>+oo ' 



as desired. 

Lemma 2.12.4. The functions ^4+ t cind Bj^^t given in Lemma \2.12.1\ satisfy 



□ 



(2.12.7) 



\A 



+,t\ 



+ \B. 



+,t\ 



1 on 



Proof. The differential equation in Lemma (12.12.21 ) implies that |^+,tp + |^+,tP is indepen- 
dent of t. Therefore, this result follows from Lemma [2. 12.31 □ 

Each of the above lemmas has analogue for the functions ^_ t and -B-,t> which we omit. The 
above lemmas together with their analogues, give the following 

Proposition 2.12.5. Let S = (a, 6; ^i, . . . , Ci, . . . , Cm) be scattering data such that 
r = ^ and are in ff^(M). For each t G M there exist unique functions A-^t — 1, B^^t, ^+,t ~ 
and B*^ ^ in H\_ (M) such that 



(2.12.8) 

and 

(2.12.9) 

(2.12.10) 

where 

(2.12.11) 



6e2*«* a* 



-B+,t 



A~,t -Bl^t 



n_-A_,t 

a 

U^-B^^t 
a 



-U^Q+,tBl, 
^^Q+,tAlt 



QUO = E 



k=l 



Furthermore, these functions satisfy 



(2.12.12) 

and 

where 
and 



dt 



q{t) 
-q{t)* 2ii 



A±M) 
B±M) 



A±^t - 1 
B±.t 



h\m.) 








as t ^ zboo, 



q{t) = -HBl,)iO' 



-q{tr = -J^{B^,t){0+). 

TT 
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2.13. Applying the discrete algorithm to continuum scattering data 

Let S = {r;(^i, . . . ,(^m',Ci, . . . , Cm) be reduced continuum scattering data witli r, G 
The main theorem tells us that there is a unique potential g : R ^ C corresponding to S. However, 
to practically compute q{t) for some t G M, one needs to somehow discretize the Marchenko 
equation 

il + U+r;U^rt)ht = -U+rt. 

Let us describe one method of doing this. 

Choose a time step A. We can replace rjA by the periodic function 

TT ^ — ' 

n=—oo 

which is an approximation to r^A in a neighborhood of zero. Let us consider f^A as a function 
of w = e^*^'^, so fjA is in H^{S^). There is a unique solution hj/^ G H^{S^) to the discretized 
Marchenko equation 

(1 + n+r*^U_rjA)~hjA = -n+fj*. 
After solving this equation, we can approximate the potential by 



(2.13.1) 



q{jA) = -HhjA){0^) « A~%a{1)- 



TT 



This, of course, greatly resembles the inverse scattering theory for hard pulses ! Let us make the 
correspondence explicit. We have 



k=l 



which implies that 

Therefore, 

TjaH = 



fjA{2nA) = f(2(n - j)A) - 2m ^ Qe^^^'^^^-")^. 



fe=i 



A 

TT 



W 



f(2(n-i)A)w;"-2Az^Cfe ^ e^i^k{j-n)A^n 

=— oo k=l n=— oo 

. —1 m —1 

- ^ f(2(n-j)A)u;"-2Ai^Cfce2^«'=^'^ ^ 

-oo 



k=l 



m 7 

W-Wk 
K=l 



where 



f{w) = — w f(2nA)u;"', 



n=— oo 



and 



Cfc = 2AiwkCk. 

So we see that replacing the reduced scattering data S by the reduced discrete scattering data 

S = {r;wi,...,Wk;ci,...,Ck) 
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is algorithmically equivalent to discretizing the Marchenko equation in the above manner. In the 
discrete algorithm, we set 

Hj = y-^ arctan \ ^j\ 

where = hjA(l)- This leads to the approximation 

q{jA) « A-Vj, 

which is very close to the right hand side of (I2.13.1I ) when 7j is small. 

In light of the above discussion, it may seem as though the discrete Marchenko equation is noth- 
ing more than a simple discretization of the continuum Marchenko equation. Such a discretization 
has been discussed in the literature, for example see [7]. However, we need to consider the following 
subtlety, which explains why the discrete theory is needed to obtain good results. Discretizing the 
left and right continuum Marchenko equations, separately, in the above sense, will not produce the 
correct pulse. Instead, one should first replace the scattering data by discrete scattering data, in the 
above manner, and then derive the data for the left equation. This will guarantee that the resulting 
hard pulse has the correct scattering data. In particular, the reflection coefficient corresponding to 
the discrete potential will be a very good approximation to the original reflection coefficient in a 
neighborhood of zero. 
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CHAPTER 3 



Pulses with finite rephasing time and applications 



In most NMR applications the designed pulses have a fixed rephasing time, p. In this chapter 
we derive a simple, SLR-type algorithm for generating hard pulses with finitely many rephasing 
time steps. We then describe several applications in NMR pulse design. 

3.1. A recursive algorithm for pulses of finite rephasing time 

Part (b) of Theorem 11.6.11 tells us that designing a hard pulse with a fixed number of rephasing 
time steps, p, amounts to specifying a function, ro, which is meromorphic in the unit disk and 
vanishing at the origin. Once r = w^frQ has been specified, one can, of course, generate the pulse 
using the recursion described in Section 12.41 However, there is a simpler, more direct recursive 
algorithm which can be used in the case of finite rephasing time. This algorithm, which we derive 
below, resembles the SLR algorithm (see [12]). 

We are deaUng with a potential of the form 

lit) = E H^it-j^)- 

j=-oo 



It is easy to check that 



So let us set 



ro = iim = — — -. 

j^oo A- J A-^p 



Notice that R^ j is a meromorphic function on the unit disk which vanishes at the origin. The 
recursion (12.8.91 ) induces the following recursion on R-j: 

^ 1 - -/jW-^R-j+i 
7* + w-^R-j+i ■ 

Since R^ j vanishes at the origin, we must have 

R-,j+i 



7j 



w 



w=0 



Thus we can reconstruct the potential from the initial data -R-,p = tq. 
For example, we can specify ro as the ratio of two polynomials: 

r - R - ^P^""^ 
with Pp{w) = 0. Then, the recursion is simply 

where7. = |jf}. 
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R- 



This very much resembles the SLR recursion. Notice, however, that there is no need to choose 
polynomials A and B satisfying |j4p + |i?p = lon the unit circle. As a result, the resulting pulses 
will generally have infinite duration. 

3.2. Equiripple pulse design 

3.2.1. The 1ST method. In this section, we describe a method for designing a pulse with a 
fixed rephasing time p < oo, which gives a profile which uniformly approximates some real ideal 
magnetization profile. Suppose we are given an ideal profile 



MideaK^:) 



_2ndcal 
1+lndealP 



l+kidealP 



where ndeai : M ^ M is a real reflection coefficient. According to Theorem 11.5. li b), we should 
uniformly approximate rjdeai by a real reflection coefficient r = e~^*^''ro where tq has a meromor- 
phic extension to the upper half plane and lim|^|_^oQ ro(^) = 0. In fact, we should design tq to 
be analytic in the upper half plane, so that the resulting pulse has minimum energy (see Corollary 
11.7.21) . So the Fourier transform of r should be supported on the half ray [— 2p, oo). Since r is to be 
real, J^{r) should actually be supported on the symmetric interval [—2p, 2p]. Therefore the problem 
reduces to uniformly approximating a real function rjdeai by a function whose Fourier transform is 
supported on a given interval [—2p, 2p]. To practically implement this procedure it is best to work 
in the discrete theory, and use an algorithm such as the Remez algorithm. 

Let us focus on the most typical example where ndeai = which corresponds to a single 

shce selective 90° pulse. To use the Remez algorithm, the user would specify the time step A, and 
three of the following parameters: 

(i) The rephasing time: p; 

(ii) The transition width: r; 

(iii) The in-slice ripple: (^ijst; 

(iv) The out-of-slice ripple: (52,ist- 

The unspecified of these four parameters can then be determined using, for example, the parameter 
relations given in 113|. The Remez algorithm then produces a periodic function (of period ^) 
which approximates rideai with a maximum error (^ijsT inside the interval [—1, 1] and a maximum 
error of 52,ist outside of the interval [— 1 — r, 1 + r]. The algorithm does not attempt to control 
the function in the transition region [— 1 — r, — 1] U [1, 1 + r]. See Section [3.2.3| for plots of pulses 
obtained using this method. 

3.2.2. The SLR method. The SLR method is a procedure for designing a pulse with duration 
T such that the resulting flip angle profile approximates some ideal flip angle profile. The duration 
is controlled by specifying B from part (c) of Theorem 11.5.11 The Fourier transform of B must be 
supported on [0, 2T] C M, where T is the desired pulse duration. This function is designed so that 
1 — 2|i?p approximates the z-component of the desired magnetization profile (or the cosine of 
the flip angle profile). One can then compute A to be analytic and non-vanishing in the upper half 
plane with |Ap = l — |i?ponM. The reflection coefficient is given by 



r(0 



A{0 

By specifying the rephasing time, p, and the zeros of A in the upper half plane, one has some limited 
control on the phase of the transverse magnetization. 

Let us focus on the case of a selective 90° pulse where ndeai = X[-i,i] and |i?ideai| = ^X[-i,i]- 
Again, for practical purposes, it is best to work in the discrete theory so that B is a. polynomial. To 
design this polynomial, the Remez algorithm can be used with the following parameters: 
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(i) The rephasing time: p; 

(ii) The transition width: r; 

(iii) The in-slice ripple: 5i,slr; 

(iv) The out-of-slice ripple: 52,slr- 

As before, three of these parameters are specified by the user, and the fourth is determined by the 
parameter relations from 1,13 J . 

3.2.3. Comparison of the SLR and 1ST methods for selective 90° pulses. In this section we 
compare the SLR and 1ST methods from Sections 13.2.11 and 13.2.21 We make the comparison for 
various values of the rephasing time, p, the transition width, r, and the out-of-slice ripple, (^2,trans- 
Here 52,trans represents the maximum magnitude of the transverse magnetization, M^; + i'M.y, for 
out-of-slice frequencies (recall that M^; + i'M.y is ideally zero out-of-slice). One can check that 
'52,trans is related to (52,ist and (52,slr by 

2(^2 1ST 

(3.2.1) 02,trans = ^ ^2 



(3.2.2) <52,trans = 2(52,SLR y 1 - (^f^sLR- 

Given r, p, and 52,trans, the in-slice ripple (5i,iong is determined. This is defined to be the maximum 
error in the longitudinal magnetization, M^, for in-slice frequencies (recall that is ideally zero 
in-slice). One can check that ^i^ong is related to 5i,ist and (^i,slr by 



(3.2.3) (5i,iong 



1 - (5l,IST - ^<Ji,ist) 
2 



(3.2.4) (^i,iong = 2^/25i,slr + 25? SLR- 

In Figures 13.2.11 13.2.21 and 13.2.31 we compare SLR and 1ST pulses for various values of the 
transition width, r, the rephasing time, p, and the out-of-slice ripple 52,trans- We see that, in each 
case, the inverse scattering pulse produces a better profile. Of course, the 1ST pulses are somewhat 
longer in duration. In many applications, however, the extra duration causes no problem, because 
the important duration is often the duration of the portion of the pulse following the peak (see IIIOI ). 

3.3. Self refocused pulse design 

Suppose we want to design a pulse with zero rephasing time (p = 0). Such a pulse is called a self 
refocused pulse. According to part (c) of Theorem [T3TT] we should approximate the ideal reflection 
coefficient, rideai, by a function r : M ^ C which has a meromorphic extension to the upper half 
plane with limi^i^oo r{z) = 0. For example, r could be a rational function with numerator degree 
strictly smaller than denominator degree. Of course, there are many ways to approximate rjdeai in 
this way. The energy formula (11.7.11 ) tells us that the energy of the resulting pulse depends on the 
locations of the poles of r in the upper half plane. If energy is a major concern, then r should be 
designed to have a small number of poles, close to the real axis. In this way, there is a delicate trade- 
off between the energy of the pulse and the accuracy of the approximation. Another concern is the 
stability of the pulse under imperfect magnetic field conditions. For example, in many applications 
it is desirable for the pulse to maintain its selectivity when it is scaled by, say, 90% or 1 10%. 

In this section we describe one method for designing relatively low energy, self refocused 90° 
pulses. We need to approximate rjdeai = X[-2,2] by a function r which has a meromorphic extension 
to the upper half plane. Let us write 

^ = 

where i? : M ^ C has an analytic extension to the upper half plane. The idea is to choose R, so that 
|e^| is very large in-slice and very small out of slice. Then r will be close to rijeai- The magnitude 
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Figure 3.2.1. Comparison of SLR and 1ST equiripple pulses: p = 3, (52,tia 
0.01, r = 0.2 X 27r. 
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of is determined by the real part of R, so we should design Re i? to be a smooth function of the 
form 

if|e|<2 
^-^2 if|e|>2 + 2r, 

where ki, k2, and r are positive numbers. The imaginary part of R is then chosen so that R has an 
analytic extension to the upper half plane. The parameters ki and k2 control the in-slice and out-of- 
slice errors, respectively, and r is the transition width. One can experiment with different values of 
these parameters to obtain a variety of pulses, with different energies. We plot one of these pulses 
in figure 13.3.11 

Remark: Notice that the r designed in the previous paragraph does not satisfy lim|^|_^oo r(^) = 
0. This can be easily remedied by subtracting an appropriate, very small constant. 



3.4. Half pulse design 

In some applications it is only necessary to approximate the x-component of the ideal magneti- 
zation profile using a self refocused pulse. For example, see lUTl . Let Mjdeai.x : M ^ M be an ideal 
x-magnetization profile. We want to find a function r, which has a meromorphic extension to the 
upper half plane with lim|^|^oo r{z) = such that 

2Rer 



1 + \r\ 



M, 



The trick is to set 



1-R 
l + R' 
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Figure 3.2.2. Comparison of SLR and 1ST equiripple pulses p 
0.1, r = 0.2 X 27r. 
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Figure 3.2.3. Comparison of SLR and 1ST equiripple pulses: p = 1, (52,trans 
0.05, r = 0.5 X 27r. 
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where 



Then we have 



R= . 

1 + r 



2Re r r + r* 



1 + \rr 1 + r*r 



l-R I 1-R* 
_ 1+R 1+R* 

~ 1 _i_ 1-R 

^ l+R* 1+R 

(1 -R)il + R*) + (1 + - R* 
(1 + i?)(l + -R*) + (1 - R){1 - R* 

1 + \R\^' 

Therefore, we should design R so that 

where M^; is a perhaps smoothed out version of M.idsai,x- The above calculations lead to the follow- 
ing 

Proposition 3.4.1. Let M.^ : M ^ M Z^e an x -magnetization profile with sufficient smoothness 
and decay, and assume that \'M.x{z)\ < 1 for all z G M. Then there exist infinitely many self 
refocused pulses u; : R ^ C such that Tuj = M. These pulses are parameterized by the poles of 

in the upper half plane, where r is the reflection coefficient. 

Remark 3.4.2. The special case of no poles is often a minimum energy pulse. For example, 
this is the case whenever M^; is non-negative. Indeed then we have \R\ < 1, which implies, by the 
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Figure 3.4. 1. A half pulse and its resulting magnetization 
2i 1 i ; 1 1 




fisqusncy (2 pi) 

maximum modulus principle, that R is never equal to — 1 in the upper half plane, which is equivalent 
to saying that r = has no poles in the upper half plane. 
Figure [3AT] shows a pulse designed using this method. 
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CHAPTER 4 



Conclusion 



We have seen that the hard pulse approximation (the approximation of an RF-pulse by a sum 
of (5-functions) leads to a discrete scattering theory which is completely analogous to the standard 
continuum scattering theory for the ZS-system. We introduced the DIST algorithm, a recursive 
algorithm for solving the full discrete inverse scattering problem relating to NMR pulse design, and 
we explained how this algorithm could be used to efficiently approximate the continuum inverse 
scattering transform. In the past, numerical techniques have been used to approximate the solutions 
to the Marchenko equations. In this thesis we have provided a more exact method of pulse design, 
which involves replacing the standard scattering data by, so called, discrete scattering data. 

The case of pulses with finite rephasing time is particularly simple, and useful in practice. We 
described some new applications to NMR pulse design. Specifically we explained how to maintain 
control on the phase of the magnetization profile while producing equiripple pulses, self-refocused 
pulses, and half pulses. 
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APPENDIX A 



Notation and background material 

A.O.l. Functions on the real line and the unit circle. In this thesis, we work with complex 
functions on the real line, M, and the unit circle, S*^ . If / : M ^ C is a function on the real line, then 
we say that / has an analytic (meromorphic) extension to the upper half plane, H, if there exists an 
analytic (meromorphic) function / on BI such that 

lim f{x + iy) = f{x) 

for almost every x G M. Such an extension is necessarily unique. In this case we abuse notation 
and let / also denote this extension /. Similarly, if / : 5^ ^ C is a function on the unit circle, then 
we say that / has an analytic (meromorphic) extension to the unit disk, D, if there exists an analytic 
(meromorphic) function / on D such that the radial limits exist almost everywhere, and coincide 
almost everywhere with /. Again, we let / denote both the function on and its extension to the 
unit disk. 

When they are not being used to represent complex numbers, the symbols ^ and w will denote 
the identity functions on M and S^, respectively. For example, if r : M C is a function on the real 
Une, then so is ^r. Furthermore, if r has a meromorphic extension to the upper half plane, then so 
does ^r. 

Let a* denote the complex conjugate of the complex number a. If / is a function on the real 
Une, we let /* denote the complex conjugate of /. If / has an analytic (meromorphic) extension to 
the upper half plane, then we consider /* also as an analytic (meromorphic) function on the lower 
half plane given by 

no = ficr- 

We follow similar conventions if / is a function on S^. Specifically, if / has an analytic (meromor- 
phic) extension to D, then /* is considered as an analytic (meromorphic) function on C \ D given 
by 

Here, C is the Riemann sphere. 

A.0.2. The Fourier transform. For > 0, we let L'^'(M) denote the space of equivalence 
classes of measurable functions / : M ^ C for which the quantity 

/oo 
-oo 

is finite. The functions / and g in L^(R) are considered to be equivalent if ||/ — ^'H^/b =0. The 
space (M) is a Hilbert space with respect to the inner product 

/oo 
f{z)g*{z)dz. 
-oo 
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We let T : L^(IR) L^(M) denote the Fourier transform, and we write ^(/) = /. For integrable 
/, the Fourier transform takes the form 



/oo 
/( 
-oo 



The Fourier transform is a unitary map, meaning it is a linear isomorphism which preserves the 
inner product, up to a factor of 27r. We let J^~^ denote the inverse of J^. For integrable / we have 



1 r°° 

f{z)=J^-\f ){z) = - j J{t)e^'^dt. 



Similarly, we let L^{S^) denote the space of measurable functions / : S*^ — > C for which the 
quantity 

II/IIl^ ■■= \f{e''td9f^ 

is finite. The functions / and g in L^{S^) are considered to be equivalent if ||/ — gW^k = 0. The 
space L'^{S^) is a Hilbert space with respect to the inner product 

We let ^ : L'^{S^) denote the discrete Fourier transform, and we write J-'{f) = f. Here 

L^{S^) is the Hilbert space of sequences / : Z — >^ C for which 



■■= E \f(j)\' 

i=-oo 

is finite. The inner product on L^{S^) is given by 

oo 

'J ^9) := E fmur 

Explicitly, we have 

/(j) = Hf){j) ■■= ^fjf{e^%e-^^'de. 

The discrete Fourier transform is a unitary map, meaning it is a linear isomorphism which preserves 
the inner product, up to a factor of 27r. Again, we let denote its inverse. For absolutely 
summable sequences S = {ajj^g^ we have 

00 

J^-\S){e^<^) = ^ a,^^'. 

j=-oo 

A.0.3. The Sobolev Space H''. We say that / : M ^ C has a weak derivative /' : M ^ C if 
for all smooth functions : M — > C with compact support we have 



f 

J —t 



f {x)(j) {x)dx + I f' {x)(p{x)dx = 0. 



The weak derivative of a function on is defined similarly. 

Let A = M or S"^, and let be a non-negative integer. The Sobolev space H'^{K) C ^^(A) 
consists of those functions which have k weak derivatives in L^. The spaces iJ^(M) and H''{S^) 
are Hilbert spaces with respect to the inner products 

^ poo 

(/'5):=E / d^fi^)d'9*{x)dx 
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and 



rr'a Jo 



j=0 

respectively. Notice that = L'^{Q). 

These spaces are particularly simple in the Fourier domain. Let us define 

H''{R) := |n G L'^{R) : H {1 + \x\f''\u{x)\'^dx < oo 



{u,v) := I + \- x^^)u{x)v* {x)dx 



with the inner product 

(iL. n) : = 
and 

( oo 

H^{S^) := i n : Z ^ C : + I^D^'^I^HI^ < °° 



with the inner product 



:= (l + n2 + ...+ 



ii?^)u(n)v{n) 



FactA.0.3. Let K = M.or S^, and let f be in {A). Then < 
Proof. By the Cauchy-Schwartz inequality we have 



1/(^)1^ 



ny)f{y)dy 



oo 

X 



< 2 / \ny)\\f{y)\dy 

J — oo 

< 2{f,f) 

< 2||/||z.2||/'|L, 

< ii/iii^ + 11/11^2 

= WfWm- 

□ 

Fact A.0.4. Let A = 'Ror S^. If f and g are in H^{A), then the product fg is also in H^{A), 

and 

||/ff|lHi< 2 11/11^1 Mh^. 

Proof. We have 

\\f9\\m = \\fg\\l2 + Ufgy\\l2 = \\f9\\l2 + \\f'g + fg'\\l, 

< ll/llio.||<7|li2+2||/'5||i. + 2||/5'||', 

< ll/llio.||<7|li2+2||/'||^,|b||ioc+2||/||i^||5'||i. 

< 2||/||i^(||5||i. + ||<7'||'.) + 2||/'||^,Nli.. 

The last inequality uses Fact |A.0.3l □ 

Fact A. 0.5. Let A = RorS^. The Fourier transform is a unitary map from H^{A) onto 
H^{A). 
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A.0.4. Self-adjoint operators. A bounded operator A : H —>^ H on a. complex Hilbert space 
is called self-adjoint if 

{Av,w) = {v,Aw) \/v,w £ H. 
It is easy to verify tiiat for sucli an operator {Av, v) is real for every v e H. We say that A is positive 
if 

{Av,v) > 

for all nonzero v ^ H. 

The following fact can be found in ll8l . 

Fact A.0.6. If A : H ^ H is a bounded, self-adjoint operator on a complex Hilbert space, 
then its operator norm is given by 

\\A\\ = sup I {Av^ v) \. 

\v\=l 

The adjoint A* of a bounded operator A: H ^ H onsi complex Hilbert space is defined by 

{Av^w) = {v,A*w) yv,w G H. 
The following two facts are easy to verify. 

Fact A. 0.7. If A and A* are bounded operators on a complex Hilbert space, then the operator 
A* A is self-adjoint and positive. 



Fact A. 0.8. Let A : X ^ X be an operator on a Banach space. If \\A\\ = a < 1, then 1 + A 
is invertible. 



Lemma A. 0.9. Suppose that A : H ^ H is bounded, self-adjoint and positive. Then 

1 
2 



A-1 



- 2' 



whenever 7 > 

Proof. This follows from Fact lA.0.6l and the estimate 

(^>(^-|)^') = {v,Av) -^{v,v) <^{v,v) 



which holds since 



< {v, Av) < \\A\\ {v, w) < 7 {v, v) . 



□ 



Proposition A. 0.10. Suppose that A : H ^ H is a bounded, positive self-adjoint operator 
on a complex Hilbert space. Then \ -\- Ais invertible, and 

||(l + ^)-i|| < 1. 

Proof. The following method was used by Epstein in |i5J. The trick is to set 7 = \\A\\ and 
write 

2 + 7, 



where 



Lemma |A.0.9| tells us that 



l + A = -^{l + B), 

2 7 

B = [A--). 

2 + 7^ 2' 

, 27 7 

LB < = — — < 1, 

' "-2+7 2 2+7 
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and so 1 + S is invertible. Clearly, 1 + A is also invertible. 
Given v e H,v/e know that 



and hence 



\{l + B)v\>{l-^)\v\, 



\{1 + A)v\ = ^l\il+B)v\ 



> ^(1 ^ 



2 ^ 2 + 7^ 

which proves that 



[1 + A)-H\ < 1. 
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APPENDIX B 



The error from softening a pulse 



As mentioned in the introduction, it is necessary, in practice, to replace a given hard pulse 



j=-oo 



by a softened version 



^^soft(t)= X ^X[o,A)(i-jA). 



J=-0O 



In this section we estimate the difference between the magnetization profiles Mhard Mgoft re- 
sulting from Ohard and Q^ofu respectively. 

Fix a frequency z G M. Let Af_,hard(^; ^j) and M__soft(^; ^j) denote the magnetizations at 
time t = Aj (or the time step). These are normalized by 



lim M_,hard(^; Aj) = lim M_,,^n{z; Aj) 

J— »— oo j-^— oo 






1 



We are interested in the difference between M_ hard(-2; jA) and M_,soft(-z; J A) for large j. 

Let us focus on the error introduced at the j* time step. Without loss of generality, we can as- 
sume that is real and positive. For the hard pulse, M_^hard(^; j + is obtained from M_^hard(^^; j) 
by a rotation of Uj radians around the y-axis, followed by a rotation of Az radians around the z-axis. 
On the other hand, one can check that M_ soft(^; j + 1) is obtained from M_ soft('Z;i) by a single 
^ ' " " 

rotation of 7 



OJ? + (Az)^ radians around the 



Wj/7 
Az/j 



-axis. The difference between these 



rotations is a rotation of the sphere (around some axis) of some number 9j of radians. If we can 
estimate Oj, then we can estimate the maximum spherical error introduced at the time step. 
Let Ry be rotation by loj radians around the y-axis, let Rz be rotation by Az radians around the 



z-axis, and let Ryz be rotation by 7 radians around the '^j/'f 

Az/j 



-axis, where 



7 



a;2 + (Az)2. 



The question is: How well does Ryz approximate RzRy"^ 
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These rotations are represented in SU2 by 



Rz = 



Ryz — 



COS 



sm ■ 



sm^ 

COS 



COS 4^ + i sin ^ 



cos ^ + sin ^ 



COS 



z sm 



- — Sin 77 

7 2 



cos ^ — sin ^ 



Using a symbolic mathematics computer program, we compute the real part of the upper left com- 
ponent of the matrix Ry^R^Ry G SU2 to be 



4>{u;j,Az) : 



.UJj. .Az. ,7, LUj UJ.j Az, ,7, 

cos(y ) cos( — ) cos(-) + — sm(y ) cos( — ) sm(-) 

Az ,ujj. . ,Az. . ,7, 
+— cos(^) sm(— ) sm(-), 



which implies that Ry^R^Ry is a rotation of 2 cos ^ (t){ojj,Az) radians around some axis. From 
the Taylor expansion 



1 



1 



1 



(cos-^ 4>{a,P)Y = -u:%Azf - ^^H^^f - ^^|(^^)' + • • • ' 

we see that, for reasonably small uoj and Az, Ryz is extremely close to RzRy composed with a 
rotation of radians. In fact, numerical evidence suggests that 2cos^^ (j){u}j,Az) < 4 



2 I'^J 
for all oj and Az 



UliZ 



Remark B.0.11. The above conclusion would be the same if we replaced the axis of rotation 
of Ry by some other axis orthogonal to the z-axis. 

Applying the above result, the maximum error between Mhard(^) and Msoft(2;) introduced at 
the j^'^ step is conjectured to be a rotation of at most ^I'^i^l radians, which would indicate a total 
error of at most 



E = 



A\z\ 



J=-0O 



AUI 



\^softit)\dt 



radians on the sphere. Thus we see that by reducing the time step. A, the error introduced by 
softening a pulse can be made as small as we like over a fixed frequency interval. 
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APPENDIX C 



The case of non-simple zeros 



In this section we consider the case where the poles of a in the upper half plane (or unit disk) 
are not necessarily simple. We simply outline the idea, and leave most of the details to the interested 
reader. 

In the case of simple zeros, the bound state data for the continuum scattering transform is 
encapsulated in the rational function 

m „ 



For each t G M, the rational functions 



Q-AO = 



e-6 



k=l 

can be uniquely determined from (5+,o using the following properties: 



(1) Q+,t has the form Q+,t(0 = E 

(2) has the form Q_,t{0 = Ek=i ^t^' 



(3) Q+,t = n_Q+,,e2^«(*-^) for all t > s; 

(4) Q_^t = n_Q_,se-2^^(*-*) for all t < s; 

(5) The function aQ-^tQ+,t + ^ is analytic in the upper half plane. 

The fifth property can be verified by observing that the residue of aQ-^tQ+,t + ^ Cfe is 

CkCka'i^k) + = 0. 

« (Cfc) 

Let us now handle the case where a has non-simple zeros. Let ^i, . . . , be complex numbers 
in the upper half plane, and let di, . . . , be positive integers (representing the multiphcities). In 
this case, we specify the bound state data by defining the rational function 

where Pk is a polynomial of degree dfe — 1 which does not vanish at ^k- For each i G M, the rational 
functions Q+^t and Q-^t are defined by the following properties: 

(1) Q+^t has the form (5+,t(C) = J2^=i ^^-^f ' where Pk,t is a polynomial of degree A; — 1; 

(2) Q-^t has the form Q~,t{0 — YlT=i ^^'^^^^^^ ^k,t is a polynomial of degree A; — 1; 

(3) Q+^t = n_Q+,,e2^€(*-*) for all t > s; 

(4) qj^t = n_Q_',,e-2i?(t-^) for all t < s; 

(5) The function aQ-^tQ+,t + a i^ analytic in the upper half plane. 
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Once these functions have been determined, one can define the right and left Marchenko equations 
using 

St = n_se-2^«*-Q_,t. 
A similar method can be used in the discrete case. 
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APPENDIX D 



Explicit implementation of DIST 



Let us write down some explicit formulas so that the DIST recursion can easily be implemented 
on a computer. We leave the derivations to the reader. 
Given discrete scattering data 

S = {a,b;wi, . . . ,Wm;c[, . . . ,c'^), 

we set 

b 



and 



r = 



6* 
a 



Cfe 



a'{wk) 
a'{wk) 



• Step 1: Define the following sequences: 

m 

f{n) = f{n) - ^CkW^'' 



k=l 

m 



g{n) = s(n) -^Cfcu;^" ^ 
Step 2: Choose M+ >> and set 





-,M+ 




' 1 " 


. ^4 


-,M+ 








Then define the polynomials and L+^j (for j = 0, 1, . . . , M+) recursively using 



jW w 



and 



-77 



^ E^=o/(-i-n)^+j(n) 
• Step 3: Choose M_ >> and set 







' 1 " 










Then define the polynomials K^j and L_ j (for j = — M_, ...,—1,0) recursively using 





Ki+1 











1 



7i 
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and 

IlZ=o9U -n)K-,j{n) 



7j 



• Step 4: Set 

M+ 
j=-M- 

where 

7* 

Uj = 2Z|-^ arctan \ 

Remark D.0.12. If the left and right values of 70 are inconsistent, then M+ and M_ should 
be increased. Initially, these integers should be chosen so that 

00 

n=M+ 

and 

00 

i-= Yl 

n=M_ 

are small. 
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APPENDIX E 



Scattering on Lie groups 



E.l. The scattering transform on Lie groups 



The selective excitation transform fits into a more general framework. Let G be a Lie group 
acting on a space X, with a special point xq G X. Fix an element J € g in the Lie algebra of G 
such that the one parameter subgroup {e"'^ : 2; € M} C G is isomorphic to S^. Also, fix a subspace 
k C g. Let a;:M^kCgbea function with sufficient smoothness and decay, so that there is a 
solution V- : M X R ^ X to the equation 

(E.1.1) dtv-{z;t) = {Lo{t) +zJ)v-{z;t) 

normalized by 

(E.1.2) lim e~-^'=*v_(z;t) 



Xq. 



For now, we assume that the solution v- is unique, and that 

{TG,.J,\<i,X,xo^){z) 



lim e'-^''*v-{z;t) G X 



exists for all z. The operator 'Tg, J,k,x,xo is called the scattering transform. It maps k- valued func- 
tions of time into X-valued functions of frequency. 





' ■ 







1 


" 




If we take G = SO^{M), X = = 





,J = 


-1 








G so3(M), and 




1 



















Im X 






-Re X 



— Im X 
Re X 




X G C 



then Tq^j^\^^x,xq coincides with the selective excitation transform T defined above. For other exam- 
ples of scattering transforms, see Section (IE.3I ). 



E.2. The scattering equation 

Let us consider the scattering transform in the special case where X 
acts on itself by left multiplication. To simplify notation we write 

^G,J,k := '?G,J,k,G,l- 
This scattering transform is universal in the following sense: 

Proposition E.2.1. Let G, J, k, X and xq be as in Section Wl\ Then 

'TG,J,k,X,xo^ = i'TG,J,kOj){xo)- 

We omit the proof, which is a straightforward application of the definitions. 

In this universal case, we consider two solutions tj-i- : M x M ^ G to the equation 

(E.2.1) dtv±{z;t) = {uj{t) + zJ)v±{z;t) 

normaUzed by 

(E.2.2) lim e--^^*w±(z;t) = 1 G G. 

4— » it 00 



G, xo = 1 G G, and G 
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The function 

(E.2.3) S^z) ■.= v+iz;t)-\^{z;t) 

is independent of t because given any two solutions vi and V2 to equation (IE.2.11 ) we have 

dtV^^V2 = -V^^{dtVl)v^^V2+V^^dtV2 

= -V^^ {uj{t) + Zj) V2 + V^^ {uj{t) + Zj) V2 

= 0. 



This calculation also demonstrates the uniqueness of v- and v^. Letting t tend to oo, we find that 

S^{z) = lim e-^^V(z;t) = {Tg,j^^uj){z) 

so that 5*1^ is the scattering transformation Tcj^k applied to u. By Proposition IE.2. II we have 

whenever G acts on a space X. 
It is useful to define 

u±{z; t) := v±{z; t)e~"'''*. 
In terms of this notation equations (IE.2.11 ) through (IE.2.31 ) become 
(E.2.4) dtu±{z; t) = uj{t)u±{z; t) + z {Ju±{z; t) — u±{z; t)J) 



lim u±{z; t) = 1 £ G, 

t— >itoo 



(E.2.5) 
and 

(E.2.6) S^{z) = e-J'*v+{z-t)-^u.{z;t)e^'K 

Equation (IE.2. 61 ) corresponds to the scattering equation (12.3.61 ). 

E.3. Three scattering transforms 

Let us fix X = C to be the Riemann Sphere with = and let e"^^ = 5^ be the group of 
rotations around the origin. There are three basic examples of scattering transforms (see Section 
EB: 

• Euclidean: G = Ge = Aut{C) is the group of invertible affine transformations 

Ge = {z^ Xz + c: A, c G C, A / 0} . 

• Spherical: G = Gs — is the group of rigid rotations of the sphere. 

• Hyperbolic: G = Gr — PSL2^ is the group of conformal automorphisms of the unit 
disk. 

The second example corresponds to the selective excitation transform. Let us show that the first 
example (Euclidean) corresponds to the inverse Fourier transform. 
We can represent 



It is natural to choose k 



and 



Ge = 
c 





1 
. Let us write 



U- {z; t) 



L0{t) 





1 U-{z;t) 
1 
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Then equation (IE.2.41) is 

dt 

or 



1 u^{z;t) 
1 



uj{t) 




1 U-{z;t) 
1 



+ z 



iu^{z;t) 




dtU- {z; t) = uj{t) + izu- {z; t). 



The solution normalized at — oo is 



U- {z; t) 



and so 



Suj{z) ■ xo= uj{s)e '■''''ds. 

J — oo 

Therefore Tg, j,k,x,xo is the inverse Fourier transform. 

E.4. The discrete scattering transform on Lie groups 

Let G, J, X, and xq & Xh& as in Section lEU That is, G is a Lie group acting on X, and J G g 
is an element of the Lie algebra such that the one-parameter subgroup e'''^ C G is isomorphic to 
S^. This time we consider functions of the form Q, : Z ^ K C G, where K is some subset of 
G. Let us identify with the subgroup e""^ C G, and set Qj = Suppose there is a solution 

V- : X Z ^ X to the recursion 



normalized by 



V-iw;j + 1) =wnjV.{w;j) 
lim 'w~^V-{w; j) = xq. 



Again, we assume that the solution is unique, and that 

exists for every w. The transform T^'j^ x xo called the discrete scattering transform. It maps 
sequences in K to loops in X. 













1 









If we take G = SOsiR), X = R^, xq = 





,J = 


-1 








G S03(l 


), and 




1 


















K = exp 






Im X 






-Re X 



— Im X 
Re X 




X G C 



then Tq^j ^ X xo coincides with the discrete selective excitation transform T'*'**'^ defined above. No- 
tice that K consists of the rotations around axes orthogonal to the z-axis. For other examples of 
discrete scattering transforms, see Section (IE.6I ). 

E.5. The discrete scattering equation 

Let us consider the discrete scattering transform in the special case where X = G, xq = 1 is 
the identity element, and G acts on X by left multiplication. To simplify notation we write 

/T-disc /T-disc 

^G,J,K ■— ^G,J,K,G,l- 

This discrete scattering transform is universal in the following sense: 
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Proposition E.5.1. Let G, J, K, X and xq be as in Section WA Then 

We omit the proof, which is a straightforward application of the definitions. 

In this universal case, we consider two solutions V± : x ^ G to the recursion 

(E.5.1) V±{w-j + 1) = wnjV±{w-j) 

normalized by 

(E.5.2) lim w-^V±{w]j) = 1 E G. 

j— >±oo 

The function 

(E.5.3) Si''{w):=V+{w-3r^V^{w;3) 

is independent of j because given any two solutions Vi and V2 to the recursion ( IE.5.11 ). we have 

Vi{w;j + l)-^V2{w-j + l) = {Vi{w;j)-^n-^^w-^){w^,V2{w-j)) 

= Vi{w;j)-^V2{w;j). 

This calculation, in particular, shows that V- and Vj^ are unique. Letting j tend to 00, we find that 

Sl{^\w) = lim w~W^{w-j) = {Ti^l^m^) 

SO that is the discrete scattering transformation T^'j^ applied to Vt. 
It is useful to define 

U±{w;j) := V±{w;j)w~^. 
In terms of this notation, equations (IE.5.11 ) through ( IE.5.31 ) become 

(E.5.4) U±{w;j + 1) = wnjU±{w; j)w~'^ 

(E.5.5) lim U±{w;j) = 1 G G, 

i->±oo 

and 

(E.5.6) S^^'iw) = w-W+{w-j)-^U-{w-j)wK 

Equation (IE.5.61 ) corresponds to the scattering equation (I2.4.7I ). 

E.6. Three discrete scattering transforms 

As in Section IE.3I , we fix X = C to be the Riemann Sphere with xq = and let e"^^ = 
be the group of rotations around the origin. There are three basic examples of discrete scattering 
transforms (see Section |R4l ): 

• Euclidean: G = Ge = Aut{C) is the group of invertible affine transformations 

Ge = {z^ Xz + c: A, c G C, A / 0} . 

• Spherical: G = Gs — SO^R is the group of rigid rotations of the sphere. 

• Hyperbolic: G = Gh — PSL2R is the group of conformal automorphisms of the unit 
disk. 

The second example corresponds to the discrete selective excitation transform. Let us show that the 
first example (Euclidean) corresponds to the discrete inverse Fourier transform. 
We can represent 



Ge = !^ q ^ : AgC*,xgc|. 
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It is natural to choose K 



and 



Then equation (IE.5.41) is 



or 



1 C 
1 



. Let us write 

1 UJi 



U-{w;j) 



1 



1 U-{w-j) 
1 



1 U.{w;j + 1) ■ 




' w ' 




" 1 UJj' 




1 




1 




1 





1 U^{w;j) 
1 



U-{w;j + 1) = w{U.{w;j) +ujj). 
The solution normalized at — oo is 

U-{w;j) = ^ UkW' 

fc=— oo 

and so 



A;=— oo 



Therefore Tq^j ^ ^.^ is the discrete inverse Fourier transform. 
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